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Abstract 


This  report  examines  the  electromechanical  effects  in  ice.  This  group  of  physical 
phenomena,  which  was  found  and  studied  relatively  recently,  broadens  basic 
knowledge  of  ice  and  may  have  some  practical  applications.  The  electrome¬ 
chanical  phenomena  in  this  monograph  are  separated  into  three  groups:  1 ) 
Effects  in  which  electromagnetic  fields  are  generated  by  means  of  mechanical 
actions  such  as  elastic  stress,  plastic  strain,  fracture  or  friction;  2)  Effects  in 
which  an  application  of  electric  fields  modifies  such  mechanical  properties  of 
ice  as  its  plasticity,  elasticity  and  friction;  3)  Effects  in  which  plastic  strain 
changes  electrical  conductivity  or  dielectric  permittivity  of  ice.  Experimental 
results  and  theoretical  models  are  discussed  and  some  possible  practical 
applications  suggested. 
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FOREWORD 


At  the  present  time,  thousands  and  thousands  of 
people  around  the  world  deal  with  ice,  snow  and  perma> 
frost.  They  are  scientists,  educators,  engineers,  naviga¬ 
tors,  meteorologists  and  others.  While  a  small  fraction 
of  these  people  contribute  to  the  knowledge  base  in  ice 
physics,  all  of  them  use  knowledge  from  it  frequently. 
Moreover,  successful  applied  research  is  based  upon 
fundamental  science — one  more  reason  for  ice  special¬ 
ists  to  have  a  textbook  on  ice  physics  on  their  desks. 

The  first  modern  ice  physics  text  was  Fletcher’s 
book  on  The  Chemical  Physics  of  Ice  (1970).  Fletcher’s 
book  is  in  typical  textbook  format:  it  is  reasonably  brief 
and  easy  to  understand.  He  touched  on  a  few  of  the 
most  important  topics,  but  not  all  of  them. 

The  most  comprehensive  book  on  ice  physics  to 
date  was  published  by  Hobbs  in  1974.  Hobbs  consid¬ 
ered  almost  all  of  the  basic  aspects  of  ice  as  understood 
at  that  time.  Moreover,  he  described  and  compared  sev¬ 
eral  (sometimes  opposing)  viewpoints.  This  fundamen¬ 
tal  and  rather  large  (837  pages)  book  is  commonly 
known  as  the  ‘Tee  Bible”  by  specialists  in  the  field.  In 
1974  and  1975,  two  CRREL  Monographs  on  ice  were 
produced  by  John  Glen.  These  were  briefly  and  clearly 
written  and  reviewed  almost  all  ice-related  subjects. 
This  work  was  (and  in  some  respects  still  is)  a  magnifi¬ 
cent  introduction  to  ice. 

Finally,  in  1981  Maeno  wrote  a  simple,  popular 
book  for  the  express  purpose  of  attracting  people’s  at¬ 
tention  to  the  subject. 

During  the  past  20  years,  a  significant  amount  of 
new  experimental  and  theoretical  work  has  appeared, 
dramatically  changing  our  views  on  ice  physics.  As  a 
result,  we  are  now  able  to  formulate  physical  laws  using 
more  simple  and  direct  methods.  We  have  found  some 
of  the  physical  models  used  in  the  past  to  be  completely 
wrong.  The  physics  of  ice  is  a  much  better  developed 
subject  than  it  was  20  years  ago. 

For  the  above  reasons,  I  feel  the  time  is  ripe  for  a 
contemporary  book  on  ice  physics,  incorporating  the 
known  and  proven  with  almost  20  years  worth  of  mate¬ 
rial  not  covered  by  previous  works. 

I  have  tried  to  prepare  a  “readable”  book,  and  not 
one  that  requires  the  reader  to  be  a  uniquely  educated 


person.  It  is  my  intent  to  present  the  material  in  such  a 
way  that  any  reader  attracted  by  the  title  Ice  Physics 
will  be  able  to  comprehend  it.  This  is  quite  difficult  for 
a  book  dedicated,  not  to  a  particular  field  of  knowledge, 
but  to  a  specific  material.  Indeed,  for  ice  it  means  we 
have  to  consider  a  wide  variety  of  subjects,  including 
quantum  chemistry,  solid  state  physics,  the  theory  of 
elasticity,  ionic  conductivity,  synchrotron  x-ray  topog¬ 
raphy,  crystal  growth,  the  physics  of  surfaces  and  more. 

The  primary  goal  is  to  produce  as  simple  a  book  as 
possible  without  sacrificing  scientific  accuracy.  Experi¬ 
mental  facts,  physical  ideas  and  theories  will  be 
strongly  organized  and  bound  together  cohesively.  The 
reader  will  be  introduced  to  a  wide  variety  of  material 
on  a  step  by  step  basis.  Then  the  picture  will  be  whole. 

To  accelerate  materials  publication,  this  book  will 
appear  first  in  the  form  of  a  series  of  joint  CRREL- 
Dartmouth  reports,  later  to  be  published  in  CRREL’s 
Monograph  series,  on: 

1.  The  structure  of  ordinary  ice 

Part  I:  “Ideal”  structure  of  ice.  Ice  crystal  lattice 
Part  II:  Defects  in  Ice 
Volume  1 :  Point  defects 
Volume  2:  Dislocations  and  planar  defects 

2.  Electrical  properties  of  ice 

Part  I:  Conductivity  and  dielectric  permittivity  of 
ice 

Part  II:  Advanced  topics  and  new  physical  phe¬ 
nomena 

3.  Optical  properties 

4.  Electro-optical  effects  in  ice 

5.  Thermal  properties 

6.  Mechanical  properties  of  ice.  Elasticity  and  anelas- 
tic  relaxation.  Plastic  properties.  Fracture  of  ice 

7.  Electromechanical  phenomena  in  ice 

8.  Surface  of  ice 

9.  Other  forms  of  ice  and  their  properties 

10.  Ice  in  space 

11.  Ice  research  laboratories 

The  reports  will  be  prepared  in  a  sequence  conve¬ 
nient  to  the  author.  The  present  report  is  the  sixth  in  the 
series. 
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Electromechanical  Phenomena  in  Ice 


VICTOR  F.  PETRENKO 


INTRODUCTION 

Among  the  numerous  interesting  properties  of  ice, 
there  is  a  group  of  physical  phenomena  in  which  elec¬ 
trical  and  mechanical  properties  correlate.  That  corre¬ 
lation  can  manifest  itself  in  three  ways: 

1.  Electromagnetic  fields  are  generated  by  means 
of  mechanical  actions,  such  as  elastic  stress,  plastic 
strain,  fracture  or  friction. 

2.  Mechanical  properties,  such  as  plasticity,  elas¬ 
ticity  and  friction,  are  modified  by  the  application  of 
electrical  fields. 

3.  Electrical  conductivity  or  dielectric  permittivity 
are  changed  by  strain. 

We  will  refer  to  such  effects  as  electromechanical 
phenomena  (EMP),  Most  of  the  experimental  results 
described  in  this  report  were  discovered  recently,  dur¬ 
ing  the  last  decade.  Of  course,  ice  is  not  unique  in 
exhibiting  electromechanical  effects.  Similar  phenom¬ 
ena  were  previously  found  in  semiconductors,  ferro¬ 
electric  and  piezoelectric  materials,  where  EMP  may 
have  a  much  greater  intensity — electrical  fields  are 
much  stronger,  for  instance.  Nevertheless,  EMP  in  nat¬ 
ural  ice  may  have  quite  a  large  total  magnitude,  such  as 
electrical  potential  difference,  owing  to  the  large  size 
typical  of  ice  sheets.  In  addition,  since  ice  is  one  of  the 
most  widespread  materials  on  Earth,  the  electrome¬ 
chanical  effects  described  below  may  have  important 
practical  applications. 

ELECTRICAL  PHENOMENA 
IN  ICE  FRICTION 

This  section  describes  two  recently  discovered 
effects  in  friction  on  ice  and  snow.  The  first  one  is  a 
strong  frictional  electrification  (Petrenko  and  Colbeck 
1995),  and  the  second  is  an  effect  of  electrical  fields  on 


ice  friction  (Petrenko  1994a).  These  phenomena  illus¬ 
trate  the  important  role  that  ice  electrical  properties 
plays  in  ice  friction.  It  is  most  likely  that  these  effects 
result  from  the  special  structure  and  electrical  proper¬ 
ties  of  the  ice  surface.  Since  the  physical  properties  of 
the  ice  surface  were  discussed  in  detail  in  a  previous 
report  (Petrenko  1994b),  we  will  only  briefly  review 
the  information  necessary  to  the  following  discussion. 

Structure  and  electrical  properties 
of  the  ice  surface 

The  physical  properties  of  the  ice  surface  differ  sig¬ 
nificantly  from  those  of  bulk  ice.  In  the  temperature 
range  from  to  approximately  -4°C,  there  is  a  thin 
liquid  or  liquid-like  film  on  the  ice  surface,  which  has 
optical  properties  very  similar  to  those  of  ordinary 
water  (Furukawa  et  al,  1987,  Elbaum  et  al.  1993).  The 
film’s  thickness  depends  on  temperature,  surface  crys¬ 
tallographic  orientation  and  ice  purity.  The  film  is  also 
sensitive  to  environmental  conditions  and  it  is  about  50 
nm  at  ~1°C.  A  special  layer  remains  on  the  ice  surface 
even  at  temperatures  down  to  -100°C.  This  layer  exhi¬ 
bits  a  diffusion  coefficient,  electrical  conductivity, 
structure  and  viscosity  that  are  quite  distinct  from 
those  of  bulk  ice  (see  Maeno  1973,  Barer  et  al.  1977, 
Golecki  and  Jaccard  1978,  Mizuno  and  Hanafuza 
1987). 

One  of  the  best-known  theoretical  models  devel¬ 
oped  to  explain  ice  surface  structure  and  properties 
was  invented  by  Fletcher  (1962,  1963,  1968).  The  key 
idea  in  Fletcher’s  model  is  the  assumption  that  at  the 
ice  surface  a  significant  fraction  of  molecules  (denoted 
as  tto)  is  oriented  with  their  dipole  moments  (i.e.,  pro¬ 
tons)  pointing  outwards,  as  opposed  to  the  bulk  where 
molecules  are  oriented  randomly.  This  would  result  in 
the  buildup  of  a  large  positive  polarization  charge  at 
the  surface.  Fletcher  came  to  his  conclusion  about 


molecule  orientational  ordering  by  analyzing  the  inter¬ 
action  between  dipole  and  quadruple  moments  of  wa¬ 
ter  molecules  close  to  the  surface.  Such  an  orientation 
of  water  molecules  implies  a  high  density  surface  elec¬ 
trical  charge  Though  this  charge  was  predicted  by 
Fletcher  many  years  ago,  its  existence  was  proven 
experimentally  only  recently  during  experiments  on 
the  frictional  electrification  of  ice  and  snow. 

All  other  predictions  of  Fletcher’s  model  are  conse¬ 
quences  of  the  molecular  ordering  in  the  surface  layer. 
In  fact,  since  in  the  bulk  there  is  no  ordering  between 
the  surface  and  the  bulk,  there  should  be  a  transitional 
layer  with  a  large  concentration  of  hydrogen  bond  de¬ 
fects  allowing  molecules  to  reorient.  This  transitional 
disordered  layer,  according  to  Fletcher,  is  the  liquid¬ 
like  layer.  Besides  the  Fletcher  model,  there  have  been 
several  others  designed  to  explain  the  unique  physical 
properties  of  the  ice  surface.  No  one  of  them  is  able  to 
account  for  all  of  its  known  properties. 

One  important  surface  property  is  ice  friction.  The 
bibliography  on  ice  and  snow  friction  is  vast  (see  re¬ 
views  by  Hobbs  [1974]  and  Colbeck  [1992]).  Here,  we 
can  examine  only  the  major  ideas  available  in  the  liter¬ 
ature  on  this  subject.  In  doing  so,  we  will  underscore 
the  facts  that  make  ice  different  from  most  other  mate¬ 
rials.  In  static  friction,  the  strong  and  universal  adhe¬ 
sion  of  ice  to  almost  any  solid  is  one  such  difference. 
Depending  on  temperature,  ice/slider  interface  rough¬ 
ness  and  slider  materials,  it  takes  from  a  fraction  of  a 
second  to  hours  to  produce  mechanically  a  very  strong 
ice/solid  interface.  Intensive  mass  transport  via  the  liq¬ 
uid-like  surface  film  accelerates  this  process  at  temper¬ 
atures  above  -10°C.  The  nature  of  the  strong  bonding 
between  ice  and  solids  is  not  quite  understood  yet  and 
perhaps  originates  from  the  special  arrangement  of 
water  molecules  on  the  ice  surface  mentioned  above. 

The  coefficient  of  kinetic  friction  of  ice  in  the  tem¬ 
perature  interval  from  -3  to  -40°C  strongly  depends 
on  sliding  velocity  v,  varying  from  0.9  at  v  =  10“^  m/s 
to  0.05  at  V  >  10"^  m/s  (Barnes  et  al.  1971,  Evans  et  al. 
1976,  Jones  et  al.  1991).  For  granite-on-ice  friction, 
Barnes  et  al.  (1971)  distinguish  three  ranges  of  sliding 
velocity  characterized  by  the  different  mechanisms 
that  govern  the  friction:  ice  creep,  when  v  <  10"^  m/s; 
ice  plastic  flow  and  fracture,  when  10"^  m/s  >  v  >  10“^ 
m/s;  and  frictional  melting,  when  v>  10“^  m/s,  T  = 
-11.75°C.  The  velocity  at  which  the  transition  to  fric¬ 
tional  melting  takes  place  may  depend  on  the  initial  ice 
temperature,  normal  pressure  and  ice  and  slider  rough¬ 
ness,  but  the  occurrence  of  melting  and  liquid  lubrica¬ 
tion  of  the  ice  friction  have  been  proven  many  times. 
Yet,  for  a  precise  quantitative  description  of  such  fric¬ 
tion,  we  need  more  data  on  both  the  film  thickness  and 
the  real  contact  area.  At  temperatures  close  to  the  melt¬ 


ing  point  and  at  high  normal  pressure,  the  pressure 
melting  of  ice  can  also  contribute  to  the  formation  of  a 
liquid  film  between  the  ice  and  a  slider. 

The  study  of  the  surface  potential  of  ice  and  electri¬ 
cal  charges  at  the  ice  surface  has  quite  a  long  history 
and  continues  to  attract  the  attention  of  numerous 
scholars.  Such  interest  in  the  electrostatic  properties  of 
the  surface  of  ice  is  stimulated  by  both  the  fundamen¬ 
tal  problems  associated  with  it  and  certain  practical  is¬ 
sues.  Fundamental  questions,  which  the  study  of  the 
surface  potential  and  ice  surface  might  help  to  resolve, 
include  the  issue  of  the  structure  of  the  surface  of  ice. 
For  instance,  if  all  the  molecules  at  the  ice  surface  are 
oriented  “proton-outwards”  as  Fletcher’s  model  sug¬ 
gests,  this  must  lead  to  a  positive  charge  at  the  surface 
X^  and  a  positive  potential  of  the  ice  surface  (p^.  Both  X^ 
and  cps  are  determined  by  the  thickness  of  the  surface 
layer  A  and  the  fraction  of  the  oriented  water  mole¬ 
cules  contained  in  this  layer.  Thus,  measurement  ofX^ 
and  (ps  in  conjunction  with  other  data  (A  for  example) 
might  provide  valuable  information  on  the  micro¬ 
scopic  structure  and  nature  of  the  ice  surface  layer. 

The  practical  issues  stimulating  the  study  of  the 
electrical  properties  of  the  ice  surface,  include,  first  of 
all,  the  problem  of  atmospheric  electricity  (see  review 
in  Hobbs  1974)  as  well  as  the  problem  of  adhesion  and 
friction  of  ice  and  snow  (Petrenko  1994a,b). 

Among  physical  mechanisms  suggested  to  explain 
ice  surface  charge  and  thunderstorm  electricity  were 
surface  electronic  states  (Buser  and  Jaccard  1978),  a 
thermoelectrical  effect  (Latham  and  Mason  1961, 
Latham  1963),  charge  separation  on  freezing  (Work¬ 
man  and  Reynolds  1949,  1950),  and  the  motion  of 
charged  dislocations  produced  by  rubbing  one  piece  of 
ice  against  another  (Takahashi  1969a). 

Early  works  on  asymmetrical  rubbing 

Until  recently,  asymmetrical  rubbing  of  ice  on  ice 
was  the  only  phenomenon  studied  in  connection  with 
ice  frictional  electrification.  Chalmers  (1952)  mea¬ 
sured  an  electrical  charge  from  small  snow  fragments 
when  two  handfuls  of  snow  were  rubbed  together, 
finding  a  negative  electrical  charge  of  2-8  x  10“^^  C  on 
the  ice  fragments.  Yoshida  (1944),  Reynolds  et  al. 
(1957),  Latham  and  Mason  (1961)  and  Latham  (1963) 
studied  charge  separation  during  asymmetrical  rub¬ 
bing  of  ice  on  ice.  In  such  an  experiment,  the  constant 
rubbing  point  becomes  warmer  than  the  variable  rub¬ 
bing  point.  The  electrical  potential  differences  gener¬ 
ated  in  such  experiments  were  very  small  and  did  not 
exceed  tens  of  millivolts.  Such  weak  electrification  is 
attributable  to  the  rubbing  piece  being  heated  more 
than  the  rubbed  one,  with  the  potential  difference  aris¬ 
ing  from  the  thermoelectrical  effect.  Since  the  thermo- 
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electrical  power  of  ice  does  not  exceed  4  mV/°C  (Bry¬ 
ant  and  Fletcher  1965),  this  effect  is  always  small. 

Takahashi  (1969a)  measured  the  electrical  potential 
difference  between  a  vibrating  metal  plate  and  the  sur¬ 
face  of  pure  single  crystals  of  ice,  both  before  and  after 
the  ice  was  rubbed  with  another  single  crystal  or  with 
a  plane.  He  found  changes  in  the  potential  difference 
ranging  up  to  0.2-f).3  V  after  the  rubbing.  Unfortunately, 
he  did  not  realize  that  he  simply  produced  an  ice-metal 
contact  potential  or,  in  other  words,  he  measured  a  dif¬ 
ference  in  electron  work  functions  between  the  metal 
plate  and  the  ice.  The  observed  relaxation  time  was 
very  long  (hours  at-10°C)  and  had  nothing  in  common 
with  the  dielectric  relaxation  time  in  ice,  which  is  5  x 
10”^  seconds  at  that  temperature.  The  effect  that  Taka¬ 
hashi  found  is  likely  ascribable  to  a  difference  in  the 
structure  and  thickness  of  the  surface  charge  double 
layer  on  fresh  and  aged  ice  surfaces.  The  layer  slowly 
changed  by  adsorption  of  impurities  from  the  air  and 
diffusion  of  the  impurities  in  ice.  We  will  discuss  such 
an  electric  double  layer  in  the  next  section. 

Frictional  electrification 

This  writer  (Petrenko  and  Colbeck  1995)  studied 
frictional  electrification  on  cylindrical  samples  of 
polycrystalline  ice  grown  from  very  pure,  deionized 
and  degassed  water.  A  polycrystalline  ice  layer,  of 
about  2  cm  thickness  and  a  typical  grain  size  of  5-12 
mm,  was  frozen  onto  the  outside  of  a  stainless  steel 
cylinder  (Fig.  1).  The  cylinder  was  mounted  on  a  lathe 
located  inside  a  coldroom  where  temperature  could  be 


Figure  1.  Ice  friction  and  ice  electrification  meas¬ 
urements  in  the  laboratory  (top)  and  on  snow  in 
the  field  (bottom)  (after  Petrenko  and  Colbeck 
J994). 


controlled  within  ±0.2‘^C.  All  ice  samples  were  ma¬ 
chined  with  a  cutter  to  an  outer  diameter  of  10  cm  and 
then  subsequently  polished  with  sandpaper.  The  ice 
static  conductivity  was  10~^  m“^  at  7  =  -30°C 

and  10“^  at  7  =  — 10°C.  The  measurements  of 

as  were  made  both  with  the  apparatus  at  rest  and  dur¬ 
ing  experiments.  No  difference  was  found  among 
values  measured  before,  during  and  after  the  friction 
experiments.  This  showed  that  possible  local  melting 
and  refreezing  of  the  ice  surface  does  not  significantly 
change  the  ice  conductivity.  When  the  cylinder  rotated, 
a  frictional  force  changed  the  belt  tension  (Ti  -  Ti). 
Then,  the  friction  coefficient  p  could  be  calculated  as 


—  arcth 
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where  a  is  an  angle  shown  in  Figure  1. 

A  schematic  of  the  experimental  setup  is  shown  in 
Figure  2.  Reverse  rotation  of  the  cylinder  allowed  us  to 
measure  T2  instead  of  7]  Polyethylene,  aluminum  and 
stainless  steel  belts  of  2.5-cm  width  were  used  as  the 
sliders,  and  measurements  were  made  at  temperatures 
from  -5  to  -35°C  and  at  sliding  velocities  from  0.5  to  8 
m/s. 

To  determine  whether  friction  causes  melting,  a 
thin  thermocouple  (0.1  mm)  was  attached  to  the  outer 
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Figure  2.  Experimental  setup  (after  Petrenko 
and  Colbeck  1994). 
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Figure  3.  Potential  difference  V  across  75-\m-thiclc 
polyethylene  film  and  its  tension  T  (after  Petrenko 
and  Colbeck  1994).  Sliding  velocity  was  1  m/s.  The 
arrows  indicate  when  the  lathe  was  switched  on  and 
off.  Temperature  was  —19  C. 
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Figure  4.  Potential  difference  V  between  ice 
and  a  stainless  foil.  Sliding  velocity  was  1  m/s 
(after  Petrenko  and  Colbeck  1 994).  Tempera¬ 
ture  is  -30°C.  The  arrows  indicate  when  the 
lathe  was  switched  on  and  off. 
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Figure  5.  Dependence  of  the  potential  differ¬ 
ence  between  ice  and  stainless  steel  slider  on 
sliding  velocity  (after  Petrenko  and  Colbeck 
1994).  T  =  -10°C. 


side  of  the  thin  (75-[im)  metal  and  dielectric  belts  slid¬ 
ing  on  the  ice.  The  maximum  heating  by  friction  did  not 
exceed  1°C;  hence,  general  melting  of  the  ice  cylinder 
surface  did  not  take  place,  although  there  may  have  been 
some  local  melting  at  a  few  points  of  direct  ice-slider 
contact. 

Significant  frictional  electrification  was  found  for  all 
sliding  materials  used  (e.g..  Fig.  3  and  4).  The  potential 
differences  V  observed  were  comparable  both  when  the 
belt  slid  on  its  own  circular  track  on  the  ice  cylinder’s 
surface  and  when  it  slid  continuously  on  a  fresh  surface 
(spiral  track).  The  electrification  of  the  ice  surface  was 
not  homogeneous.  Normally,  it  was  possible  to  find  two 
positions  along  the  30-cm-long  cylinder  where  V  dif¬ 
fered  by  a  factor  of  two.  Since  the  belt  covered  many  ice 
grains,  contributions  from  individual  grains  were  not 
significant  and  added  only  a  small  oscillating  compo¬ 
nent  to  V. 

Typically,  the  potential  differences  generated  be¬ 
tween  ice  and  metal  were  about  300  V  (ice/Al,  -31.5  C, 
2  m/s).  Under  the  same  experimental  conditions,  the 
electrical  field  generated  by  the  friction  of  ice  on  poly¬ 
ethylene  reached  £  =  2.1  x  10®  V/m.  This  magnitude 
was  calculated  as  the  ratio  of  the  potential  difference 
across  the  film  to  the  film  thickness  L.  When  the  electri¬ 
fication  was  measured  in  the  temperature  interval  from 
-4.5  to  -35°C  and  for  sliding  velocities  from  0.5  to  8 
m/s,  the  sliders  always  received  a  positive  charge.  The 
electrification  increased  with  decreasing  temperature 
and  increasing  sliding  velocity  (Fig.  5).  The  maximum 
F  =  1.6  kV  was  found  at  the  lowest  temperature 
r  =  -35°C  and  the  highest  sliding  velocity  V  =8  m/s. 

From  the  measurements  of  a  current  /  passing 
through  the  electrometer  when  it  was  used  as  a  current 
meter,  we  found  that  the  coupled  ice-slider  acted  as  a 
charge  generator  (Fig.  6).  Moreover,  at  temperatures 
above  approximately  -12°C,  the  current  is  proportional 


-10  12 


ln[v(nn/s)] 

Figure  6.  Dependence  of  the  electric  current 
though  ice/metal  slider  interface  on  sliding 
velocity  at  three  different  temperatures  (after 
Petrenko  and  Colbeck  1994). 


to  sliding  velocity.  Thus 


/  =  ^  =  A.— «v  (2) 

dt  dt 

where  Q  -  electrical  charge 

'k  =  surface  charge  density  that  the  slider 
picked  up  from  the  ice  surface 
s  =  area  swept 
V  =  sliding  velocity 
t  =  time. 

When  the  input  impedance  of  a  measuring  device  is 
much  larger  than  the  resistance  of  ice  /?,  the  charge 
generated  flows  through  the  ice,  resulting  in  a  potential 
difference  V 


V=IR. 


(3) 


Since  R  for  ice  increases  exponentially  with  de¬ 
creasing  temperature  (see  the  Electrical  Properties  of 
Ice  [Petrenko  1993b]),  the  observed  electrification  ris¬ 
es  rapidly  at  low  temperatures.  The  charge  density 
X  found  at  T  =  -10°C  was  about  1.6  xlO“^  C/m^  and 
10"^  C/m^  at  -35°C.  An  even  higher  charge  density 
of  about  10^  C/rrP'  accumulated  at  the  ice/dielectric 
slider  interface.  Since  it  is  unlikely  that  a  slider  picks 
up  all  the  surface  charge  from  the  ice,  these  values  of  X 
appear  to  be  a  lower  limit. 

At  temperatures  below  “10°C,  the  density  of  charge 
collected  from  the  surface  becomes  nonlinearly  depen¬ 
dent  on  the  sliding  velocity  (see  Fig.  6).  While  at  T  = 
-10°C,  I  is  proportional  to  v,  at  I  increases  as 

and  at  r  =  -25°C  it  increases  as  v^. 

In  the  field  experiments,  we  used  an  alpine  ski  as  a 
slider  on  snow  (Fig.  1).  The  ski  itself  acted  as  a  flat 
capacitor  Cq  to  measure  an  electrical  field  of  strength  E 
generated  by  charge  density  A<sici  that  accumulated  on 
the  bottom  ski  surface.  A  larger  capacitor  C  =  2  pF, 
connected  in  parallel  with  Cq,  reduced  the  potential 
difference  V  to  a  magnitude  that  could  be  measured  by 
a  portable  data  logger  carried  by  the  skier  in  a  back¬ 
pack,  It  is  easy  to  show  that 


2£q£  C  +  Cq 


(4) 


where  £oand  £are  the  dielectric  permittivities  of  a  vac¬ 
uum  and  of  the  material  between  the  metal  plates  in  the 
ski  respectively.  The  reduction  factor  Cq/{C  +  Cq)  was 
5.1  X  10^. 

The  field  measurements  of  snow  frictional  electrifi¬ 
cation  were  made  in  deep  powder  snow  in  British 
Columbia.  Figure  7  shows  a  typical  record  of  V  taken 
at  -A^C  at  three  different  velocities.  The  general  char¬ 
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Figure  7.  Records  of  the  potential  difference  V 
across  a  capacitor  built  into  a  ski  at  three  slid¬ 
ing  velocities  (after  Petrenko  and  Colbeck  1994), 
Temperature  was  Skier  started  att  =  0  and 
slowed  at  the  moments  marked  with  the  arrows. 


acter  of  the  electrification  observed  in  the  field  experi¬ 
ments  agrees  well  with  the  laboratory  results.  One  of 
the  most  noticeable  differences  is  that,  in  the  field,  V 
first  becomes  negative  and  only  then  positive.  At  the 
highest  speed,  V  reached  2.2  V,  which  corresponds  to 
X=  3.6 X 10“^  C/m^.  This  charge  density  was  compara¬ 
ble  with  a  value  of  1 .6  x  10“^  C/m^  found  on  ice.  In  the 
absence  of  the  buffer  capacitor  C,  V  would  have  been 
4.4  X  10^  V  if  there  had  been  no  electrical  breakdown 
of  the  air  or  the  ski  dielectric  material.  And  this  is  not 
yet  at  saturation.  According  to  the  laboratory  results, 
we  can  expect  even  higher  voltages  at  lower  tempera¬ 
tures.  A  reader  interested  in  ski  electrification  can  find 
more  details  in  Colbeck’s  (1994)  publication. 

Possible  reasons  for  electrification  by  friction,  i.e., 
the  accumulation  of  the  electrical  charges  transferred 
from  the  ice  surface  onto  the  slider,  can  be  classified 
into  two  categories: 

1 .  The  slider  sweeps  up  the  charge  already  present 
at  the  ice  surface. 

2.  Friction  somehow  produces  a  spatial  separation 
of  the  charge,  which  is  then  swept  up  by  the  slider. 

After  analyzing  our  data  and  the  known  ice  surface 
properties,  we  favor  mechanism  1.  Figure  8  shows  a 
very  general  electrical  charge  distribution  in  ice  near 
the  surface.  In  a  thin  subsurface  layer,  a  charge  layer  of 
an  extremely  high  density  X^  is  formed.  The  major  fac¬ 
tor  leading  to  its  formation  is  the  ordering  of  dipole 
moments  of  the  water  molecules  in  the  upper  mono- 
molecular  layer,  resulting  in  polarization  Pq,  Such 
ordering  was  initially  suggested  by  Weyl  (1951)  for 
water  and  was  later  used  by  Fletcher  in  his  model  of 
the  ice  surface  structure.  Other  evidence  for  the  exist¬ 
ence  of  such  an  ordered  layer  of  water  molecules  at  the 
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Figure  8.  Schematic  dependence  of  the  space  charge  density  pon  the  dis¬ 
tance  to  the  ice  surface  x  (after  Petrenko  and  Colbeck  1994).  Scales  are  arbi¬ 
trary.  d  is  the  thickness  of  liquid-like  layer,  D  is  the  thickness  of  the  ice  layer 
rubbed  off  by  a  slider  from  the  ice  surface,  and  and  Kf  ^  are  the  length  of 
screening  by  Bjerrum  defects  and  ions  correspondingly. 


ice  surface  comes  from  the  high  adsorption  coefficient 
of  ice  at  temperatures  above  -40°C  (Ocampo  and 
Klinger  1983),  and  the  results  of  computer  simulations 
carried  out  by  Bush  and  Devlin  (1991).  At  the  same 
time  the  polarization  of  water  molecules  in  the  bidk  is 
zero.  Therefore,  the  polarization  P  drops  from  Pq  at 
the  surface  (;c  =  0)  to  0  (jc  =  oo).  This  nonuniform  polar¬ 
ization  results  in  a  surface  charge  of  a  high  density 
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dx  =  Pq. 


(5) 


An  estimation  using  eq  5  gives  A<s «  2  xl0“^  C/m^. 
Other  possible  reasons  accounting  for  the  formation  of 
could  be  adsorption  of  ions  from  the  air  and  the 
presence  of  the  surface  electronic  states. 

In  thermal  equilibrium  the  surface  charge  must 
be  neutralized  by  the  opposite  screening  charge  Age, 
equal  to  the  former  in  magnitude,  which  is  distributed, 
however,  in  a  much  thicker  screening  layer.  The 
screening  of  the  surface  charge  in  ice  was  considered 
in  Petrenko  and  Ryzhkin  (1984a),  and  in  Petrenko  and 
Maeno  (1987).  The  screening  charge  is  composed  of 
protonic  charge  carriers  (Bjerrum  D-  and  L-defects 
and  H30‘^  and  OH“  ions).  There  are  two  screening 
lengths  iniceKf^  and  K2^ 
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where  Eq  =  dielectric  permittivity  of  vacuum 

Eco  =  3.2,  the  high-frequency  permittivity  of  ice 
Cg  K  100  =  static  dielectric  permittivity  of  ice 
k^  =  Boltzmann  constant 
^3  =  0.38^  =  charge  of  Bjerrum  defects 
n\Q  =  concentrations  of  H30'^  ions 
n2Q  =  concentrations  of  OH~  ions 
«30  =  concentrations  of  D-defects 
rt40  =  concentrations  of  L-defects. 

At  the  screening  length  Kf  1,  which  is  determined  by  the 
concentration  of  majority  carriers  (/130  +  ^40)  and  coin¬ 
cides  with  the  well-known  expression  for  the  Debye 
screening  length  Kf^  the  field  drops  down  from  the 
original  value  As/Eoe^  to  Ag/Eoes*  As  we  have  seen,  £s  ~ 
10^,  that  is,Kf^  is  a  characteristic  length  at  which  £s 
attains  a  steady-state  value.  At  larger  x  the  field  drops 
down  to  zero  because  of  screening  by  minority  charge 
carriers,  i.e.,  ions  with  a  larger  characteristic  screening 
length  determined  by  ion  concentrations  (nio  +  n2o)* 
Substitution  of  charge  carrier  concentrations  in  pure  ice 
(Petrenko  1993b)  into  eq  7  and  8  gives  at  7  =  -5°C 

—  =  3.4xl0"^m  (8) 


—  =  3xl0“^m.  (9) 

K2 
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But  at  -35°C  they  are  longer  than  eq  8  and  9  predict 
owing  to  the  exponential  dependencies  of  the  charge 
carrier  concentrations  on  temperature.  At  -35°C 

—  =  8xl0~^m  (10) 

—  =  1.6xl0'^m.  (11) 

K2 

Notice  that  both  the  liquid-like  layer  and  the  princi¬ 
pal  screening  length  are  very  small.  That  makes  it 
likely  that  a  slider,  when  rubbing  ice,  sweeps  off  the 
surface  charge  together  with  the  most  of  the  screening 
charge.  That  leaves  just  a  small  fraction  of  all  on  the 
ice  surface  (see  Fig.  8).  This  may  explain  why  the  den¬ 
sity  of  the  surface  charge  taken  away  by  the  sliders  (X  < 
10"^  C/m^)  was  much  smaller  than  the  theoretically 
predicted  X^~2x  10”^  C/m^.  Another  factor  decreas¬ 
ing  ?iis  that  the  real  contact  area  is  much  smaller  than 
the  slider  surface,  and  hence  the  slider  does  not  sweep 
off  the  whole  area  of  the  interface  as  assumed  in  eq  2. 

The  described  mechanism  of  electrification  by  fric¬ 
tion,  thus,  though  hypothetical,  is  based  on  firm  and 
clear  deductions  from  the  surface  physics  of  ice  and 
can  be  tested  experimentally  in  the  future,  since  it  im¬ 
plies  characteristic  dependencies  on  temperature,  slid¬ 
ing  velocity,  thickness  of  the  “erased”  layer  and  dop¬ 
ing  of  ice. 

Other  mechanisms  that  may  in  principle  contribute 
to  the  frictional  electrification  of  ice  are  the  motion  of 
charged  dislocations  during  the  plastic  deformation  of 
ice  subsurface  regions  (Takahashi  1969a,  Petrenko  and 
Whitworth  1983),  a  charge  separation  by  microcracks 
in  ice  during  cleavage  of  surface  layers  (Petrenko 
1993a)  and  the  Workman-Reynolds  (1949,  1950)  ef¬ 
fect  (charge  separation  during  refreezing  of  ice).  It  is 
unlikely  that  any  of  these  phenomena  is  able  to  account 
for  the  large  magnitude  of  the  potential  difference,  1.6 
kV,  observed  in  our  research.  Under  the  most  favorable 
conditions,  the  first  two  mechanisms  generated  V  val¬ 
ues  that  are  much  less  then  1  V  and  the  Workman-Rey¬ 
nolds  effect  generated  about  100  V.  Also,  the  Work¬ 
man-Reynolds  effect  reaches  a  maximum  in  doped  ice 
samples  while  frictional  electrification  decreases  with 
doping.  As  discussed  in  the  previous  section,  the  ther¬ 
moelectric  effect  can  add  only  a  negligible  contribu¬ 
tion. 

The  fact  that  at  higher  temperatures  (-5  to  -10°C) 
the  electrification  is  proportional  to  the  sliding  veloci¬ 
ty  (V  oc  v)  while  at  lower  temperatures  (<  -25°C)  V  is 
proportional  to  can  be  explained  in  terms  of  the 
length  of  the  electrically  charged  track  on  the  ice  sur¬ 
face  Lt  that  a  slider  leaves  behind.  When  the  slider  is 


moving  at  constant  velocity,  any  electrical  charge 
picked  by  it  has  to  relax,  flowing  over  the  distance  to 
the  location  of  the  opposite  electrical  charge  behind  the 
slider.  Then  the  resistance  R  in  eq  3  is  proportional  to 
=  Td  •  V  and  V  is  proportional  to  v^.  (Tq  is  the  Debye 
relaxation  time  of  the  ice.)  Measurements  of  Tq  re¬ 
vealed  that  Lf  =  0. 1  mm  at  -5°C  but  Lt  ==  1  cm  at  -*35°C. 
Hence,  at  the  high  temperature  the  discharge  takes  place 
locally  under  the  slider  since  is  much  less  than  the 
slider  dimensions  (2.5  cm).  In  that  case  R  does  not  de¬ 
pend  on  V  and  V  ©c  v.  But  when  at  the  low  temperatures, 
Lf  becomes  comparable  with  the  slider  dimensions,  and 
/?  oc  oc  V  and  V  oc  v^. 

Effect  of  electrical  fields  on  ice  friction 

Earlier  (Petrenko  1994a),  it  was  demonstrated  that 
the  application  of  an  external  dc  bias  to  the  ice/slider 
interface  can  double  the  force  of  dynamic  and  static 
friction  between  ice  and  metals  and  between  ice  and  di¬ 
electrics.  The  experimental  techniques,  the  ice  samples 
and  experimental  conditions  used  were  similar  to  the 
ones  described  in  the  previous  section  (see  Fig.  1  and 
2).  The  measurements  were  performed  at  temperature 
intervals  from  -5  to  -30°C,  with  sliding  velocities  from 
0.5  to  8  m/s  and  the  dc  bias  within  the  range  from  -3  to 
3  kV.  All  slider  materials  created  a  strong  increase  in 
apparent  friction  coefficient  when  the  bias  IV I  >  1  kV 
was  applied.  Owing  to  the  high  dc  impedance  of  pure 
ice,  the  electrical  power  consumed  from  the  power  sup¬ 
ply  was  very  small  (about  lO"^  W  at  -30°C  for  metal 
sliders  and  even  less  for  dielectric  ones). 

Figure  9  depicts  the  changes  in  polyethylene  belt 
tension  T2  when  a  3-kV  bias  was  applied  to  the  ice  cyl¬ 
inder  and  the  foil  electrode  attached  to  the  outer  surface 


Figure  9.  Changes  in  tension  of  the  polyethyl¬ 
ene  belt  (T2)  when  a  3-kV  bias  is  switched  on 
and  off  (after  Petrenko  1994a),  Temperature  is 
-30°C.  Sliding  velocity  is  2  m/s.  The  friction  co¬ 
efficient  p  ^0.3  atV  =  0  and  p  ^0.5  atV  =  3  kV. 
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of  the  dielectric  film,  as  shown  in  Figure  1.  In  a  static 
electrical  field,  ice  behaves  as  a  conductor  so  that  the 
entire  bias  is  applied  across  the  dielectric  belt  and  the 
air  gap  between  the  belt  and  the  ice  cylinder.  Since  T2 
increases  when  frictional  force  increases,  the  effect 
shown  in  Figure  9  corresponds  to  greater  friction  with 
the  applied  voltage.  Within  the  limits  of  experimental 
error  (+20%,  determined  by  oscillations  of  T2),  there 
was  no  difference  in  the  effects  of  positive  and  nega¬ 
tive  biases  on  the  apparent  ice  friction  coefficient  p.  At 
higher  temperatures  (from  -5  to  about  -10°C),  (T2  - 
T[)  becomes  comparable  with  oscillations  in  the  belt 
tension  because  of  the  small  magnitude  of  ice  friction. 
This  results  in  the  fact  that  while  the  effect  of  a  dc  bias  is 
noticeable  at  these  temperatures,  it  becomes  much  more 
pronounced  only  when  the  temperature  decreases. 

The  change  in  the  frictional  force  can  be  explained 
by  means  of  electrostatic  pressure  Pgl  -  which  increases 
the  normal  pressure  of  the  belt  on  ice 
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(12) 


where  V  =  applied  voltage 

E  =  electrical  field  strength  at  the  ice/slider  in¬ 
terface 

£0  _  dielectric  permittivity  of  a  vacuum 
L  =  thickness  of  the  belt. 

Pel  reaches  a  maximum  in  the  absence  of  an  air  gap 
between  the  ice  and  the  slider.  In  this  case  and  under 
conditions  indicated  in  Figure  9,  Pgl  “  ^  kPa,  which 
exceeds  the  average  normal  pressure  Pq  of  the  belt  on 
ice  by  a  factor  of  four.  Under  the  experimental  condi¬ 
tions  described,  we  can  expect  a  fivefold  increase  in 
the  apparent  friction  coefficient.  Since  we  see  a  smaller 
effect,  we  assume  that  there  was  an  average  (and  dynam¬ 
ic)  air  gap  between  the  ice  and  the  slider  of  about  0.2 
mm,  sufficient  to  reduce  the  effect  on  the  observed 


magnitude.  Also,  Pel  should  not  depend  upon  the  bias 
polarity,  in  agreement  with  the  experimental  results. 

When  metallic  belts  made  of  aluminum  or  stainless 
steel  foils  were  used,  the  effect  of  dc  bias  on  friction 
was  comparable  to  the  results  obtained  with  dielectric 
sliders  (see  Fig.  10).  In  this  case  the  dc  bias  was  ap¬ 
plied  across  the  entire  thickness  of  the  ice  between  the 
central  metal  cylinder  and  a  metal  slider.  The  effects  of 
positive  and  negative  bias  on  friction  were  similar  to 
each  other  but  not  quite  reversible.  Namely,  when  the 
bias  was  switched  off,  it  took  a  long  time  (several  min¬ 
utes  at  -30°C)  to  restore  the  initial  frictional  force. 
Since  the  electrical  field  vanishes  during  the  much 
shorter  time  of  about  0. 1  second  (this  was  a  time  con¬ 
stant  of  the  electrical  circuit  used),  the  observed  effect 
of  the  electrical  field  on  ice-metal  friction  cannot  be 
explained  in  terms  of  electrostatic  pressure.  However, 
even  more  serious  evidence  in  favor  of  disregarding 
the  importance  of  electrostatic  pressure  comes  from  a 
quantitative  comparison  of  Pgi  and  Pq.  While  P 0  re¬ 
mained  the  same,  as  in  experiments  with  dielectric 
sliders  (about  1  kPa),  Pgi  has  not  exceeded  0. 1 8  Pa  (see 
eq  12).  This  smaller  electrostatic  pressure  comes  from 
a  larger  L,  across  which  the  dc  bias  V  was  applied.  In 
this  case  L  was  equal  to  the  ice  thickness  (1.5  cm)  since 
the  central  metal  cylinder  and  the  metal  belt  acted  as 
two  electrodes  attached  to  a  poor  conductor  (ice).  The 
absence  of  any  sizable  potential  drop  at  the  ice/metal 
interfaces  was  proven  by  perfectly  linear  current  volt¬ 
age  characteristics.  The  field  effect  was  more  pro¬ 
nounced  at  low  temperatures  for  the  same  reason  as  in 
the  case  of  the  dielectric  sliders,  that  is  better  signal-to- 
noise  ratio. 

The  application  of  electrical  fields  also  increased 
static  friction  between  ice  and  metals  (see  Fig.  11). 
This  phenomenon  was  very  noticeable  at  temperatures 
above  -15°C,  at  which  it  takes  just  a  few  seconds  for  a 
metal  belt  to  freeze  to  the  ice  cylinder.  At  lower  tem¬ 
peratures,  when  a  slider  rested  just  few  seconds  be- 


Figure  10.  Changes  in  tension  of  the  aluminum  belt 
(Ti)  when  a  3-kV  bias  is  switched  on  and  off  (after 
Petrenko  1994a).  The  smaller  Ti  is,  the  larger  fric¬ 
tion  is  (see  Fig.  1).  Temperature  is  -30°C.  Sliding  vel¬ 
ocity  is  2  m/s.  The  friction  coefficient  p  ^0.34  atV  = 
0  and  p  -0.47  immediately  after  V  =  3  kV  is  turned 
on. 


8 


10 


Figure  11.  Changes  in  tension  of  the  stainless 
steel  foil  belt  when  the  ice  cylinder  direction  of 
rotation  reverses  (after  Petrenko  1994a).  Tem¬ 
perature  is  ’-10°C.  Peaks  correspond  to  static 
friction.  The  two  sets  of  dashed  lines  indicate  the 
average  magnitude  of  the  belt  tension  Ti  and  T2 
when  2-kV  dc  bias  between  ice  and  steel  was 
*'on.  **  The  static  friction  coefficient  |X  -0.267  at  V 
=  0and\i-^0.312atW  =  2kV. 


€ 

(0 


1 

f-1 

[■“I 

— 

— ■ 

--1 

L 

L 

"1 

r 

H 

H 

t— 

L_L 

V  =  2kV 


20  40  60  80  100 

Time  (s) 


tween  sliding,  there  was  usually  no  pronounced  differ¬ 
ence  between  static  and  dynamic  friction,  and  the  elec¬ 
trical  field  affected  friction  in  the  manner  shown  in 
Figure  10.  Once  again,  a  strong  field  effect  such  as  that 
shown  in  Figure  11  cannot  be  explained  in  terms  of 
electrostatic  pressure. 

As  working  hypotheses  accounting  for  the  phenom¬ 
ena  observed,  I  considered  electrostatic  pressure  (see 
eq  12),  the  effect  of  static  electrical  fields  on  the  plastic 
deformation  of  ice  and  expansion  of  the  contact  area 
caused  by  electrical  fields. 

Electrostatic  fields  externally  applied  to  the  inter¬ 
face  between  ice  and  a  dielectric  slider  are  able  to  sig¬ 
nificantly  alter  normal  pressure  and,  consequently,  the 
friction  force  (eq  12).  Perhaps  a  more  important  and 
interesting  fact  is  that  such  electrostatic  attraction  be¬ 
tween  ice  and  a  slider  always  occurs  by  itself  owing  to 
the  frictional  electrification.  Under  certain  conditions, 
such  electrification  can  contribute  to  the  formation  of 
the  friction  force. 

The  effect  of  electrical  fields  on  plastic  deformation 
of  thin  ice  monocrystals  (Petrenko  and  Schulson  1993) 
was  used  to  explain  the  action  of  an  external  electrical 
field  on  the  dynamic  friction  of  metal  sliders  against 
ice.  It  appears  that  an  application  of  relatively  small  dc 
electrical  biases,  on  the  order  of  50  V,  to  an  ice/metal 
interface  reduces  by  nearly  an  order  of  magnitude  the 
plasticity  of  a  layer  of  ice  of  50-100  pm  thickness  ad¬ 
jacent  to  the  metal.  This  in  turn  is  caused  by  the  extrac¬ 
tion  of  the  major  charge  carriers  (Bjerrum  defects) 
from  the  region  near  the  electrode  by  the  electrical 
field  (Petrenko  and  Schulson  1992).  Depletion  of  Bjer¬ 
rum  defects  from  the  subsurface  layer  of  ice  alters  the 
conditions  of  dislocation  motion,  which  results  in  the 
reduction  of  plasticity.  When  ice  slides,  there  is  a  plas¬ 
tic  shear  deformation  in  the  contact  zones  between  the 
ice  and  slider.  Therefore,  a  reduction  in  the  plasticity 
of  ice  when  a  voltage  Vis  applied  results  in  an  increase 


of  resistance  to  shear  stress,  or  in  other  words  results  in 
an  increase  of  frictional  force. 

The  effect  of  an  electrical  field  on  the  plasticity  of 
ice  may  also  contribute  to  the  observed  effect  of  a  field 
on  static  friction  (see  Fig.  1 1).  In  addition  to  that  mech¬ 
anism,  there  is  the  possibility  that  the  applied  field 
changes  the  contact  area  between  the  ice  and  metal. 
The  following  grounds  support  this  assumption:  the  ef¬ 
fect  of  an  electrical  field  on  static  friction  was  ob¬ 
served  only  when  the  temperature  was  high  enough 
that  a  quasi-liquid  film,  characterized  by  low  viscosity 
and  a  large  self-diffusion  coefficient  as  compared  to 
ice,  existed  at  the  ice  surface.  When  a  strong  electrical 
field  is  present  in  the  voids  between  ice  and  metal,  ice 
is  attracted  into  these  regions  by  the  electrostatic  pres¬ 
sure  and  tends  to  fill  them.  Since  the  static  dielectric 
permittivity  of  ice  is  quite  large  (65  ~  100),  the  driving 
force  of  electrostatic  pressure  is  also  quite  large 

(13) 

A  mass  transfer  at  high  temperatures  can  efficiently 
occur  along  the  quasi-liquid  film,  filling  air  gaps  be¬ 
tween  the  ice  and  the  metal  slider  and  increasing  the 
ice  adhesion. 

The  above  experimental  facts  and  discussion  show 
the  importance  of  electrical  effects  on  ice  friction. 
Physical  mechanisms  of  the  phenomena  are  not  deter¬ 
mined  yet.  We  have  considered  here  only  a  few  obvi¬ 
ous  mechanisms  that  are  known  to  be  applicable  under 
similar  experimental  conditions  and  can  be  easily  esti¬ 
mated.  Other  mechanisms,  such  as  electrolytic  effects 
caused  at  the  metal/ice  interface,  the  role  of  thin  local 
liquid  patches,  modification  of  the  mechanical  properties 
of  ice  surface  by  multiple  local  melt-refreeze,  etc., 
should  also  be  analyzed.  This  is  an  extensive  field  for 
future  study. 
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ELECTRO-ELASTIC  EFFECTS 

Is  ordinary  ice  I^  piezoelectric? 

Perhaps  the  most  well-known  electromechanical 
phenomenon  is  piezoelectricity.  Application  of  uni¬ 
form  elastic  strain  to  a  piezoelectric  crystal  causes  its 
electrical  polarization.  The_^electrical  polarization  is 
characterized  by  a  vector  P .  Crystalline  quartz  is  a 
well  known  example  of  a  piezoelectric  material.  Math¬ 
ematically,  the  direct  piezoelectric  effect  is  described 

Pi  =  ^oXij^j  "*■  ^ijk  ^  ^jk 

where  Pi  =  component  of  the  polarization  vector 
^ijk  “  piezoelectric  stress  tensor 
^ij  =  electrical  susceptibility  tensor 
=  strain  tensor. 

The  inverse  effect  produces  a  contribution  to  the  strain 

£mn  ”  ^kmn  ^k 

where  Cmnij  =  stiffness  tensor 
aij  =  stress  tensor 
Jkmn  =  piezoelectric  strain  tensor. 

Two  piezoelectric  tensors  relate 

^kmn  ”  ^kij 

where  ^ninij  is  the  compliance  tensor. 

The  physical  mechanism  of  piezoelectric  polariza¬ 
tion  is  very  simple  and  can  be  explained  in  the  following 
way.  The  atomic  structure  of  all  piezoelectric  crystals 
may  be  thought  of  as  consisting  of  pairs  of  positive  and 
negative  ions.  Since  one  such  pair  makes  an  electric 
dipole,  we  can  say  that  the  whole  structure  of  a  piezo¬ 
electric  crystal  is  composed  of  such  elementary 
dipoles.  Now,  using  this  simplified  interpretation  of 
piezoelectric  polarization,  we  can  easily  estimate  its 
absolute  magnitude,  which  can  later  be  compared  with 
the  magnitude  of  pseudo-piezoelectric  effects  in  ice. 
Let  us  take  a  typical  length  of  an  elementary  dipole  a  = 
10”io  and  the  dipole  concentration  N  =  3x  10^^  m. 
Then 

(17) 

£o£ 

Application  of  a  strain  Ea  varies  the  length  a  by  the 
factor  of  Eii.  The  piezoelectric  polarization  is  described 
with  the  change  in  the  polarization  vector 

AP-eiiP.  (18) 

Hence,  the  components  ^ijk  of  the  piezoelectric 


stress  tensor  in  eq  14  should  be  of  the  order  of  P 

eijk-F-2xl0l0C/m2.  (19) 

The  reader  can  check  that  eq  19  gives  the  correct  order 
of  magnitude  for  strong  piezoelectric  crystals  such  as 
segnet  salt  (4.2  x  10^^  C/m^).  In  quartz,  which  is  a  weak 
piezoelectric,  about  1.5  x  10^  C/m^, 

Figure  12  illustrates  generation  of  the  electrical 
polarization  P  under  the  action  of  uniform  compres¬ 
sion.  Looking  at  the  figure,  we  can  easily  realize  that 
the  piezoelectric  effect  is  incompatole  with  a  center  of 
symmetry.  Otherwise,  the  vector  P  would  not  know 
which  direction  to  go:  to  the  left  or  to  the  right.  Accord¬ 
ing  to  Neumann’s  principle,  the  symmetry  of  any 
physical  property  of  a  crystal  cannot  be  less  than  the 
symmetry  in  the  atomic  structure  of  that  crystal.  In  our 
case  this  means  that  crystals  with  a  center  of  symmetry 
in  the  atomic  structure  cannot  exhibit  a  piezoelectric 
effect. 


a 


Figure  12.  Generation  of  piezoelectric  polariza¬ 
tion  in  a  crystal  under  the  action  of  the  uniform, 
uniaxial  compression. 


Let  us  apply  Neumann’s  principle  to  ice.  Ice  is  a 
molecular  crystal  consisting  of  water  molecules.  While 
the  centers  of  the  water  molecules  are  arranged  into  a 
regular  wurtzite  lattice,  the  molecular  dipoles  may  have 
one  of  six  allowed  orientations  and  are  not  ordered  (see 
Hobbs  1974  and  Petrenko  1993c).  Such  a  random  ori¬ 
entation  of  the  dipoles  results  in  zero  polarization  of  ice 
as  it  is  shown  in  Figure  13b.  In  other  words,  ice  has  a 
kind  of  statistical  center  of  symmetry.  An  application  of 
uniform  strain,  as  shown  in  Figure  13c,  may  cause 
some  reorientation  of  a  small  fraction  of  the  molecules 
but  without  distortion  of  the  center  of  the  symmetry 
(Petrenko  and  Ryzhkin  1984b).  The  reorientation  of 
water  molecules  in  ice  under  strain  is  responsible  for 
the  ice  anelastic  relaxation.  We  can  conclude  now  with 
confidence  that  ordinary  ice  Ih  cannot  be  an  ordinary 
piezoelectric.  Though  some  early  papers  reported  ob¬ 
servation  of  a  piezoelectric  effect  in  ice  (see,  for  in¬ 
stance,  Rossman  1950),  it  is  likely  that  this  resulted 
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b.  Chaotic  orientatioi^  of  water  c.  Partial  ordering  of  the  water  mole- 

molecules  in  ice  before  stress.  cules  in  ice  after  stress  waj  applied. 


Figure  13.  Application  of  Neumann’s  principle  to  ice. 


from  some  of  the  pseudo-piezoelectric  effects  dis¬ 
cussed  later  in  this  report. 

Phonon-induced  polarization  of  ice 

Whalley  and  Klug  (1984)  noticed  that  while  a 
macroscopic  bulk  of  ice  does  not  have  a  dipole  mo¬ 
ment,  at  the  microscopic  scale  of  about  of  tens  of  ang¬ 
stroms,  the  ice  lattice  may  be  polarized  by  statistical 
fluctuations  in  the  small  statistical  ensembles.  That 
means  that  at  the  microscopic  scale  ice  can  behave  as  a 
piezoelectric.  As  Whalley  and  Klug  showed,  this  can 
affect  propagation  of  elastic  waves  with  a  wavelength 
comparable  to  the  size  of  the  polarized  micro-regions. 
The  same  phenomenon  should  also  result  in  additional 
absorption  of  infrared  light. 

Polarization  induced  by  nonuniform  strain 

Though  ordinary  ice  Ijj  is  not  a  conventional  piezo¬ 
electric,  under  some  special  conditions  elastic  strain 
can  induce  an  electrical  polarization  in  it.  Such  special 
conditions  are  gradients  of  strain,  temperature  or  im¬ 
purity  concentration  and  also  an  externally  applied 
electrical  field.  All  these  factors  destroy  the  statistical 
center  of  the  symmetry,  whose  existence  imposes  a 
ban  on  the  conventional  piezoelectricity  of  ice. 

Let  us  start  the  consideration  of  such  pseudo-piezo¬ 
electricity  from  an  electrical  polarization  of  ice  caused 
by  nonuniform  strain.  This  phenomenon  was  first  ob¬ 
served  by  Evtushenko  et  al.  (1984),  who  have  studied 
electrical  fields  generated  by  bending  vibrations  in 
thin  pure  single  crystals  of  ice.  Later  (Evtushenko  et  al. 
1987,  Evtushenko  and  Petrenko  1991)  such  polariza¬ 
tion  was  measured  in  a  wide  temperature  range,  from  0 
to  -80°C,  on  ice  doped  with  HE,  HCl  and  NH3  The 
wide  temperature  range  and  the  doping  were  used  to 
study  the  phenomenon  of  ice  having  different  types  of 
majority  charge  carriers.  Also,  a  theoretical  model  of 
the  pseudo-piezoelectric  effect  was  developed,  which 
is  explained  below. 

It  was  suggested  that  in  ice  under  a  nonuniform 
strain,  protonic  charge  carriers  migrate  along  the  pres¬ 


sure  gradient  in  the  same  way  as  an  air  bubble  moves 
upwards  in  water  (Fig.  14).  The  reason  for  such  migra¬ 
tion  of  the  defects  is  that  their  specific  volume  differs 
from  that  of  a  water  molecule.  As  the  result,  point  de¬ 
fects  larger  then  a  water  molecule  will  migrate  in  the 
direction  opposite  to  the  pressure  gradient,  while  those 
with  a  smaller  specific  volume  will  migrate  in  the  direc¬ 
tion  of  the  gradient.  Some  important  practical  cases  in 
which  nonuniform  strain  causes  the  protonic  defect’s 
migration  and  an  associated  electrical  polarization  are 
shown  in  Figure  15.  Also,  since  there  is  an  oscillating 
gradient  of  pressure  in  a  longitudinal  wave,  such  a  wave 
is  followed  by  a  wave  of  weak  electrical  polarization  in 
ice. 

The  formation  of  any  protonic  point  defect  results  in 
deformation  of  the  surrounding  lattice.  In  the  case  of 
Bjerrum  defects,  the  defective  hydrogen  bonds  are 
weaker  than  the  normal  ones,  and  therefore  they  will  be 
longer.  This  is  obvious  in  particular  for  a  D-defect, 
where  there  exists  an  extremely  strong  repulsion  be¬ 
tween  two  protons  on  one  hydrogen  bond.  Accordingly, 
we  can  assume  that  Bjerrum  defects  “expand”  ice, 
increasing  its  volume.  If  ice  is  under  pressure  P,  then  the 
formation  energy  of  the  i''’  defect  Egi  is  increased  by  the 


Figure  14.  Air  bubble  migra¬ 
tion  against  a  pressure  gradi¬ 
ent  in  water. 
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Figure  15.  Some  practical  cases  in  which  nonuniform  strain  generates  electric  polarization  of  ice.  The 
arrows  indicate  the  direction  of  motion  of  charge  carriers. 


amount  of  work  against  this  pressure 


E^i(P)  =  E^fQ)  +  PYai 


(20)  The  quantities  Tii  are  given  by 


(23) 


where  Yai  is  the  so-called  activation  volume.*  For  the 
defects  that  increase  volume,  Yai  is  positive,  and  their 
formation  energy  is  raised  proportionally  to  the  applied 
pressure.  This  results  in  an  exponentially  decreasing 
concentration  of  such  defects  with  increasing  pressure. 
Since  the  formation  energy  of  defects  depends  on 
the  applied  pressure  P  according  to  eq  20,  in  ^e  pres¬ 
ence  of  a  pressure  gradient,  there  is  a  force  Fj  acting 
upon  the  defects 

Fj  =  -grad(Fai)  =  "Yai  grad(P) .  (21) 

To  find  a  resulting  flux  of  the  defects  jj,  we  have  to 
add  this  force  to  the  transport  equation 

7i  =[eiF-rii<I)n)ni!ai  -Dj  grad^zj  (22) 

where  €[  =  electric  charge  of  the  defects 
E  =  electric  field  strength 
n\  =  defect  concentration 
=  defect  mobility 
D{  =  defect  diffusion  coefficient 
^  =  configuration  vector 


*  Normally,  in  solids  pressure  is  not  hydrostatic,  i.e.,  an  c^2  c^3- 

Then,  instead  of  eq  20  we  should  use 

=  £^ai(0)  +  Yai(<^  1 1  +  Q2  +  Qs)  =  ^ai(O)  +  TaiC^jj  (20a) 


or 


^ai(^j)  “  ^ai(O)  tXajEjj 


(20b) 


where  is  a  deformation-potential  constant  of  type  chaige  carrier. 


rii  =  1,-1 -1,1  for  i  =  1,2, 3, 4  (24) 


and 


y  =  3.85A:b'‘oo 

where  Too  =  2.76  A  is  the  oxygen-oxygen  distance  in 
ice  (Jaccard  1964).  To  find  the  electrical  field  strength 
caused  by  the  pressure  gradient  VP,  we  should  solve 
the  system  of  eq  22-23  under  particular  initial  and 
boundary  conditions.  Thus,  in  the  case  of  one  dominat¬ 
ing  carrier  type  and  static  VP 


YaiVP_  (26) 

e.+h^.O 

while  for  two  types  of  charge  carriers,  for  example 
H3O+  ions  and  L-defects,  and  static  pressure  gradient 
the  electrical  field  is 

E  ^  (Ya4-Yal)y^  (27) 

e 

(see  Evtushenko  et  al.  [1987]  and  Evtushenko  and  Pe¬ 
trenko  [1991]).  Electrical  polarization  of  ice  caused  by 
stable  and  growing  cracks  was  calculated  in  another 
publication  (Petrenko  1993a).  Electromagnetic  emis¬ 
sion  generated  by  such  cracks  will  be  described  in  the 
section  on  electro-fracture  effects. 

Equations  26-27  show  how  the  defect’s  activation 
volumes  can  be  determined  from  measurements  of  the 
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H  (V/m) 


pseudo-piezoelectric  effect.  Indeed,  that  method  was 
used  (Evtushenko  and  Petrenko  1991)  to  determine  Yai- 
By  varying  temperature  and  doping,  we  determined 
activation  volumes  of  all  four  types  of  protonic  defects 
in  ice.  Figure  16  shows  time  dependencies  of  an  elec¬ 
trical  bias  generated  across  a  thin  ice  sample  when  a 
step-shaped  pulse  of  bending  was  applied  to  the  sam¬ 
ple.  The  ice  was  doped  with  NH3  in  such  concentration 
that  D-defects  were  majority  charge  carriers  above 
-30°C,  but  below  -35*^0  0H“  ions  dominated.  The  ob¬ 
served  decay  of  the  signals  is  mainly  attributable  to  a 
final  input  impedance  of  the  measuring  circuit  used. 
Since  the  electrical  charges  of  OH“ions  and  D-defects 
are  opposite,  so  are  the  electrical  biases  caused  by  their 
migration. 


Other  pseudo-piezoelectric  effects 

There  are  more  trivial  pseudo-piezoelectric  effects. 
They  can  be  observed  when  ice  is  already  electrically 
polarized  for  some  reason.  For  example,  because  of  the 
thermoelectric  effect  (Brook  1958),  ice  is  always  po¬ 
larized  when  there  is  a  temperature  gradient  in  it.  The 
electrical  potential  difference  generated  is  usually 
small  since  the  thermoelectrical  power  of  ice  is  low 

-^  =  ±(2-4)mV/X  (28) 

(see  Bryant  and  Fletcher  1965,  Takahashi  1966).  More 
pronounced  electrical  polarization  of  ice  can  be  in¬ 
duced  by  nonuniform  distribution  of  impurity  ions  in 
ice.  In  this  case  AV  can  reach  hundreds  of  millivolts 
and  even  1  V  (Brook  1958,  Petrenko  1992a).  Hence, 
we  can  say  that  a  piece  of  ice  acquires  an  electrical 
polarization  if  there  are  either  temperature  or  impurity 
gradients  in  it.  Elastic  strain  ey  (even  uniform)  changes 
the  dimensions  of  the  ice  sample,  changing  its  dipole 
moment.  When  a  characteristic  frequency  of  applied 


strain  0)»cq),  where  ©d  is  a  Debye  frequency,  the 
charge  carriers  in  ice  have  no  time  to  get  redistributed 
in  the  bulk  of  the  deformed  ice  and,  in  the  first  approx¬ 
imation,  the  oscillating  part  of  the  polarization  -eijp. 
At  the  lower  frequencies  ©  =  <  ©d  we  have  to  allow 
for  the  ac  bias  arising  at  the  elastic  oscillations  being 
applied  to  the  bulk  of  ice,  whose  electrical  conductivi¬ 
ty  and  dielectric  permittivity  have  a  complicated  fre¬ 
quency  dependence  (Hobbs  1974,  Petrenko  1993b). 
Such  electrical  polarization  induced  by  elastic  strain  in 
ice  with  temperature  or  impurity  gradients,  or  both, 
was  reported  by  Evtushenko  et  al.  (1987). 

ELECTROMAGNETIC  PHENOMENA 
IN  ICE  FRACTURE 

At  present  we  know  of  several  electromagnetic 
phenomena  that  are  associated  with  the  presence  or 
appearance  of  cracks  in  ice.  They  are  electromagnetic 
emissions  (EME)  from  growing  cracks  and  stable 
cracks  under  varying  stress,  and  also  crack-induced 
changes  in  electrical  conductivity  and  apparent  dielec¬ 
tric  permittivity  of  ice. 

The  most  interesting  and  intriguing  of  these  phen¬ 
omena,  EME  from  cracks,  was  primarily  found  under 
poorly  controlled  field  conditions  when  a  low  signal- 
to-noise  ratio  left  doubts  and  uncertainty  about  the 
phenomenon’s  existence.  Later,  EME  from  cracks  in 
ice  were  reproduced  and  studied  in  a  well-controlled 
laboratory  environment.  The  laboratory  results  made  it 
possible  to  develop  physical  models  that  explained  that 
electromechanical  effect.  When  the  researchers  were 
equipped  with  these  models,  they  succeeded  in  obtain¬ 
ing  reliable  and  well-reproducible  detection  of  EME 
from  cracking  in  lake  and  sea  ice.  At  the  present  time, 
the  EMP  described  below  in  this  section  are  well  un¬ 
derstood  and  have  even  been  found  useful,  both  in  ice 
mechanics  and  remote  sensing  of  stress  in  ice.  That  is 
why  we  first  start  with  a  theoretical  consideration  of 
the  physical  processes  that  cause  EME  from  cracks  in 
ice  and  other  associated  phenomena.  Then,  we  will 
discuss  laboratory  and  field  data  and  their  interpreta¬ 
tions. 

Theory 

Cracks  splitting  p re-polarized  ice 

The  main  cause  of  EME  from  cracks  in  natural  ice  is 
the  so-called  “frozen-in”  or  “intrinsic”  electrical  field 
that  is  always  presents  in  both  sea  and  freshwater  ice. 
This  electrical  field  is  generated  by  spatial  nonunifor¬ 
mity  in  the  concentration  of  ions  dissolved  in  the  ice 
bulk  (Petrenko  1992a,  1994b).  Figure  17  shows  sche¬ 
matically  such  an  electrical  field  inside  granular  ice. 
The  electrical  field  originates  from  the  nonuniform 
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Figure  17.  Schematical  representation  of  variations 
in  ions*  concentration  n,  electric  field  strength  E  and 
electric  potential  (p  along  the  crack  surface  in  granu¬ 
lar  ice. 


Such  an  electrical  field  oscillates  from  grain  to 
grain  in  the  ice  bulk  following  grain  configurations 
and  is  very  significant  even  in  non-doped  ice  grown 
from  distilled  water,  £  =  100- 1000  V/m.  In  doped  ice 
E  can  locally  be  as  large  as  10"^  V/m  (Petrenko  1992a, 
1994b). 

When  an  ice  mass  is  split  by  a  crack  rapidly  grow¬ 
ing  in  the  direction  perpendicular  to  the  electrical  field, 
the  field  generates  two  surface  charges  of  the  opposite 
sign  on  the  opposite  surfaces  of  the  crack 

^s=eo(es-l)£x  (33) 

where  E_i  is  a  component  of  the  electrical  field  perpen¬ 
dicular  to  the  crack  surface.  Inside  the  crack  the  same 
component,  exceeds  E±  by  two  orders  of  magni¬ 
tude 


(34) 

since  85  =  10^.  Taking  £j_=  10^-10^  V/m,  we  obtain  a 
potential  difference  A(p  across  a  1-mm  wide  crack 


distribution  of  impurities  in  ice.  Namely,  when  ice 
grains  grow,  impurities  are  pushed  out  by  the  phase 
boundary  from  grain  centers  to  their  peripheries,  so 
that  more  impurities  get  trapped  near  grain  boundaries. 
This  process  is  common  for  most  materials  and  is  as- 
cribable  to  lower  solubility  of  foreign  atoms  and  ions 
in  a  solid  phase  as  compared  with  a  liquid  phase.  Mo¬ 
bile  ions  trapped  in  the  ice  bulk  then  diffuse  against  the 
gradient  of  their  concentration,  generating  electrical 
charge  transport  and  an  intrinsic  electrical  field  in  ice. 
In  thermal  equilibrium,  that  diffusion  flux  Jd  is  can¬ 
celed  with  a  drift  flux  of  the  ions  7dr 

;DR+yD  =  0  (29) 

where 


Acp-lO-lOOV.  (35) 

Such  a  large  potential  difference  can  easily  explain 
why  the  intrinsic  electrical  field  dominates  among  the 
causes  of  EME  from  cracks  in  ice. 

The  processes  of  dielectric  relaxation  causes  the 
density  of  crack  surface  charges  to  decay  with  time.  In 
the  case  of  a  growing  crack,  this  density,  thus,  decreas¬ 
es  as  the  distance  from  the  crack  tip  increases.  Hence, 
eq  33-35  are  applicable  only  in  a  crack  tip’s  vicinity. 
To  calculate  far  from  the  tip,  we  should  take  into 
account  two  processes  of  relaxation  of  electrical 
charges  in  ice  (Jaccard  1964,  Gluschenkov  and  Petren¬ 
ko  1993,  Petrenko  1993a).  In  the  first  relaxation  pro¬ 
cess,  decays  exponentially  with  time  from  its  initial 
value  A.S  to 


Jdr 

(30) 

h\ 

(31) 

The  relaxation  time  of  this  fast  process  is 


where  p, 

A 

F 

e\ 


=  ion  mobility 
=  ion  diffusion  coefficient 
=  ion  concentration 
=  ion  electric  charge. 


Equalizing  the  right-hand  sides  of  eq  30  and  31,  we 
find  the  intrinsic  electrical  field 


£  =  M!vin(ni). 


(32) 


(36) 


At  r  =  -10°C  in  pure  ice,  Xy  ~  10“^  seconds  and  is 
even  shorter  in  natural  ice.  Since  the  crack  velocity  in 
freshwater  ice  varies  in  the  range  from  10^  to  10^  m/s 
(Petrenko  1992a,  Sato  and  Wakahama  1992),  the 
length  of  a  high-density  surface  charge  train  behind  a 
crack  tip  is  about  10~^  s  x(lO^-lO^)  m/s  =  0.1-1  mm. 
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The  rest  of  the  surface  charge  relaxes  very  slowly 
with  a  relaxation  time  T2 

T2=-^  (37) 

which  is  about  0.1  second  in  pure  ice.  That  means  that 
a  crack  moving  with  10^  m/s  velocity  can  have  a  100- 
m  long  electrically  charged  tail. 

Since  the  intrinsic  electrical  field  is  an  oscillating^ 
function  of  the  crack  length  (see  Fig.  17),  so  is  the  sur¬ 
face  charge  and  the  resulting  electrical  dipole  moment 
of  the  crack.  When  a  crack  grows,  its  total  dipole  mo¬ 
ment  oscillates  and  generates  an  EME.  An  average  fre¬ 
quency  of  EME/relates  to  the  crack  velocity  Vcr  and  an 
average  grain  size  g  as 


The  constants  oCaj  are  usually  called  the  deformation 
potential  constants  and  for  a  quite  wide  variety  of  mate¬ 
rials  these  constants  are  about  10  eV. 

The  force  acting  on  charge  carriers  in  the  presence 
of  a  strain  gradient  is 

F  =  -grades/)  =  a  grad(e  .  (42) 

This  force  initiates  a  flux  of  the  charge  carriers,  i.e., 
electrical  currents,  and  as  a  consequence  an  electrical 
potential  difference  appears,  as  in  the  case  of  the  bend¬ 
ing  of  ice  samples  (Evtushenko  et  al.  1984,  1987). 

To  find  the  electrical  field  in  the  area  around  a  crack, 
we  first  find  a  vector  of  polarization  P(f ,  t) .  Let  us  con¬ 
sider  ice  with  four  types  of  charge  carriers.  Then 


/  =  (38) 

8 

This  relationship  provides  an  investigator  with  a 
simple  method  of  crack  velocity  determination  using 
measurements  of /and  g  (Petrenko  1992a,  Gluschenk- 
ov  and  Petrenko  1993).  The  dielectric  relaxation  reduc¬ 
es  Xg  and  should  be  taken  into  account.  Readers  can 
find  more  details  of  the  theory  of  such  EME  in  some 
previous  papers  (Petrenko  1992a,  Gluschenkov  and 
Petrenko  1993). 

EME  from  the  pseudo-piezoelectric  effect 

The  second  important  mechanism  resulting  in  crack- 
induced  electrical  polarization  of  ice  is  the  pseudo¬ 
piezoelectric  effect  discussed  above  (see  eq  20-27  and 
Fig.  15b).  A  theoretical  description  of  such  polariza¬ 
tion  was  developed  earlier  (Petrenko  1993a). 

When  a  crack  appears  and  expands,  it  generates 
strains  in  the  surrounding  material  that  we  define  as 
ejj(r,?) ,  which  depend  upon  coordinates  r  and  time  t. 
We  have  to  find  the  electrical  field  strength  £(r,  t)  and 
electrical  potential  (p(r ,  t)  that  are  caused  by  these  elas¬ 
tic  strains.  The  charge  carriers  have  to  move  because 
their  energy  U  depends  on  strains 


f/=f/o  +  f/i(8ij).  (39) 

In  a  linear  approximation  for  an  isotropic  material 
3 

tz-t/o-ttaj^ea  =(/o-«aj-eii.  |eii|«l  (40) 
/=1 

where  we  use  the  Einstein  summation  convention 
3 

(41) 

/=! 


where  fluxes  are  determined  by  the  transport  equa¬ 
tion 


7j  =(ej£-Tij<Dn  +  aajVea)-^-Z)jVnj  (44) 

hi 

where  Gj  =  |Lij  nj  are  the  partial  conductivities  and  Z)j 
are  the  diffusion  coefficients. 

Inside  the  bulk  material 

div£  =  -divP/8«,eo.  (45) 

Taking  into  account  that  curl  E  =  Q  (quasi-stationary 
approximation)  and  that  at  infinity  (r  =  oo)  E  =  p  =  o^ 
we  can  find 


V(p  =  -E=P/eoe^.  (46) 

If  for  a  particular  crack,  the  deformation  eji(r,r)  is 
known,  then  the  system  of  eq  43^6  is  complete.  In  the 
presence  of  dielectric/conductor  or  conductor/conduc¬ 
tor  interfaces,  we  have  to  join  the  proper  boundary  con¬ 
ditions  to  this  system  of  equations.  The  system  of  eq 
43-46  has  analytical  solutions  for  many  cases  of  practi¬ 
cal  interest  that  will  be  shown  in  the  following  sections. 

Figure  18  illustrates  three  types  of  elementary 
cracks.  Notice  that  the  mode  III  crack  does  not  produce 
an  electrical  field  since,  for  an  isotropic  material,  8ii  =  0 
(see  for  example  Hellan  1984). 

For  a  mode  I  crack 


£ii  -  Ki 


(1-v) 

E 


0 

cos  — 

_ 1 

^|2nr 


(47) 
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Figure  18.  Three  modes  of  cracks. 


where  V  is  Poisson’s  ratio,  E  is  Young’s  modulus  and  Ki 
is  an  intensity  factor 


Ki  =  ^a{°°)4Td 


(48) 


where  <7(00)18  tensile  stress  at  infinity  and  /  is  the  crack 
length.  The  numerical  coefficient  P  depends  on  bound¬ 
ary  conditions  and  is  about  1.  For  the  mode  II  crack 


=  K, 


>.  Sin  — 

(l-v) _ 2_ 

E  -Jlnr 


(49) 


=  px(°o)V7t/  (50) 

where  again  P  is  about  1  and  x(oo)  is  shear  stress  at  infin¬ 
ity  (Hellan  1984). 

For  stable  cracks  or  cracks  that  move  “slowly,”  i.e., 
cracks  for  which  the  time  of  their  motion  is  longer  than 
the  electrical  relaxation  times,  the  system  of  eq  43-46 
has  a  solution  (Petrenko  1993a) 


Figure  19  shows  force  lines  of  an  electrical  field  gen¬ 
erated  by  mode  I  and  mode  II  cracks  in  ice.  Estimates 
made  using  eq  47—51  show  that  a  typieal  electrical 
potential  difference  between  infinity  and  a  crack  in¬ 
creases  from  10-100  mV  at  the  crack’s  tail  to  about  1  V 
at  the  crack’s  tip. 

The  electrical  field  generated  by  a  rapidly  growing 
crack  can  be  calculated  as 


E  =  {t,r)  =  - 


— exp(-f/Xi) 
EEq 


(52) 


t  ctct  ■ 


where  Xi  is  the  fast  relaxation  time  from  eq  36  and  the 

type  of  charge  carriers  are  majority  charge  carriers. 

Laboratory  experiments 

Today,  we  have  an  ample  collection  of  experimen¬ 
tal  data  on  the  generation  of  EME  and  electrical  sig¬ 
nals  by  cracks  in  ice.  To  the  best  of  our  knowledge,  the 
first  laboratory^' measurements  of  electrical  signals 
from  cracks  were  performed  by  Takahashi  (1969b).  He 
measured  an  electrical  charge  generated  by  the  break¬ 
ing  of  frost  under  a  temperature  gradient.  The  frost  was 
grown  either  between  two  ice  hemispheres  or  between 
two  metal  plates.  Significant  charge  separation  was 
observed  only  under  a  temperature  gradient,  which 
generates  an  electrical  field  £j_  in  the  direction  perpen¬ 
dicular  to  the  cracks  in  the  frost.  The  strength  of  elec¬ 
trical  field  Ej_  of  about  20  V/m  at  -15°C  was  generated 
by  the  thermoelectric  effect  in  ice.  The  density  of  the 
electrical  charge  separated  by  cracks  in  frost  was  esti¬ 
mated  by  Takahashi  as  10“^  C/m.  That  value  exceeds 
by  about  one  order  of  magnitude  the  charge  density 
predicted  by  eq  33.  The  discrepancy  may  be  explained 
by  the  break  in  the  frost  occurring  at  the  narrowest 
neck,  where  the  temperature  gradient,  and  hence 
VT,  exceed  their  average  magnitudes. 

Takahashi  (1983)  expanded  his  study  of  charge  sep¬ 
aration  by  cracks  in  single-crystalline  ice  samples. 
During  those  experiments,  both  sample  cross  section 
and  the  temperature  gradient  were  well  defined.  Again, 


MODE  II 

Figure  19.  Lines  offeree  of  the  electric  field  in¬ 
duced  by  stable  mode  I  and  mode  II  cracks  in  ice 
(after  Petrenko  1993a). 
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though  there  was  large  scattering  in  the  data,  signifi¬ 
cant  charge  separation  was  observed  only  in  the  pres¬ 
ence  of  a  temperature  gradient  Vr  and,  hence,  An 
average  of  that  set  of  the  results  fits  well  to  the  predic¬ 
tion  given  by  eq  33. 

Electrical  signals  from  thermocracks  in  ice  single- 
and  polycrystals  were  studied  by  del  Pennine  et  al. 
(1974).  Their  ice  specimens  were  frozen  on  the  top  of  a 
massive  copper  substrate  that  was  also  used  as  the  first 
electrode.  The  second  electrode,  connected  to  a  lock-in 
amplifier,  was  vibrated  above  the  ice  surface  to  mea¬ 
sure  the  ice-metal  contact  potential  (Kelvin  method). 
Thermocracks  were  produced  during  intensive  cooling 
of  the  copper  substrate  because  of  the  difference  in  the 
thermal  expansion  coefficients  of  copper  and  ice.  The 
cracks  in  the  ice  (which  were  assumed  to  be  mode  I 
type)  grew  in  prismatic  planes  of  ice  single  crystals, 
starting  from  the  copper  substrate  toward  the  mobile 
electrode.  Figure  20  depicts  several  such  events.  The 
potential  jumps  associated  with  crack  nucleation 
measured  from  one  to  several  volts  and  were 
characterized  by  very  long  relaxation  time.  The  relax¬ 
ation  time  varied  from  10^  seconds  at  -20°C  to  10^ 
seconds  at  -100°C.  While  the  signals’  amplitude  can 
be  explained  in  terms  of  the  pseudo-piezoelectric  po¬ 
larization  of  ice  induced  by  cracks,  their  too-long  re¬ 
laxation  time  does  not  correspond  to  any  electrical  re¬ 
laxation  time  known  for  ice.  That  time  may  corre¬ 
spond,  as  the  authors  suggested,  to  relaxation  of  elastic 
stress  in  ice. 

Fifolt  (1990,  1991)  and  Fifolt  et  al.  (1992,  1993) 
found  and  investigated  electrical  signals  from  cracks 
in  freshwater  columnar  and  granular  ice  under  uniaxial 
compression.  Temperature  gradients  and  other  possi¬ 
ble  sources  of  external  electrical  fields  were  absent  in 
these  well-controlled  experiments.  The  measurements 
were  performed  in  a  temperature  range  from  -10  to  - 
33°C.  Electrical  signals  from  cracks  appeared  as  an 


rrc) 


Figure  20,  Jumps  of  electric  potential  V  caused  by 
nucleation  of  thermocracks  (after  del  Pennine  et  al 
1974).  Open  circles  show  V  during  freezing  at  a  con¬ 
stant  rate  of  0.5  °C/min;  solid  circles  show  V  under 
warming  at  a  rate  of  0.35  °C/min. 

electrical  potential  difference  between  two  metal  elec¬ 
trodes,  which  were  either  attached  to  an  ice  specimen 
(thin  foil  electrodes)  or  placed  at  some  distance  from 
the  specimen  surfaces  (see  Fig.  21). 

Observed  signals  ranged  in  amplitude  from  4  to  70 
mV  and  lasted  from  1.5  x  10“^  to  10“^  seconds.  Most 
commonly,  the  signals  exhibited  an  exponentially 
decaying  form  (Fig.  22).  Independent  measurement 
showed  that  the  signals’  delay  time  coincided  with  the 
slow  dielectric  relaxation  time  X2  (see  eq  37).  The  signal 
type  shown  in  Figure  22  was  typical  for  small  cracks 
that  either  split  one  ice  grain  or  were  located  on  one 
intergranular  boundary.  Oscillating,  sinusoidal  type  sig¬ 
nals  were  characteristic  for  large  cracks  that  split  many 
single-crystalline  grains  or  columns  (Fig.  23). 

The  nature  of  such  signals  was  understood  after  the 
discovery  of  the  frozen-in,  intrinsic  electrical  field  in 
poly  crystalline  ice  (Petrenko  1992a).  Cracks  in  such  ice 
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Figure  21.  Signal  recovery  system  used  to  detect  emissions  (after  Fifolt  et  al  1993). 
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Figure  22.  Electrical  signal  of  a  microcrack  in  fresh¬ 
water  columnar  ice  {after  Fifolt  et  al.  1993).  T  = 
-33°Q  <5  =  6.3  MPa. 
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Figure  23.  Electrical  signal  of  a  large  crack 
splitting  many  columns  in  freshwater  colum¬ 
nar  ice  (after  Fifolt  et  al.  1993).  T  =  -30°C, 
0  =  3.78  MPa. 


split  electrically  pre-polarized  ice,  causing  the  appear¬ 
ance  of  growing  or  oscillating  electrical  dipoles  (see 
the  theoretical  section  above).  This  mechanism  was 
later  proven  in  experiments  in  which  a  strong  external 
electrical  field  was  applied  to  ice  samples  with  grow¬ 
ing  cracks  (Petrenko  1993d).  It  was  shown  that  in  this 
case  the  magnitude  of  electrical  signals  from  cracks  is 
proportional  to  the  external  electrical  field  £j_,  in  ac¬ 
cordance  with  eq  33.  An  example  of  a  signal  generated 
by  a  microcrack  in  the  presence  of  an  external  electri¬ 
cal  field  is  shown  in  Figure  24.  Again,  the  decay  time 
of  the  signals  coincided  with  the  slow  dielectric  relax¬ 
ation  time  12* 

Electrical  polarization  of  ice,  induced  by  a  stable 
crack,  attributable  to  the  pseudo-piezoelectric  effect, 
was  studied  using  the  scheme  depicted  in  Figure  25  to 
determine  electrical  field  distribution  around  the  crack 
(Petrenko  1992a).  Narrow  cracks  prepared  by  a  thin 


Figure  24.  Typical  electrical  signal  cap¬ 
tured  from  a  microcrack  in  freshwater 
ice  under  the  action  of  an  external  elec¬ 
tric  field  E_i=  10  V/cm.  T  =  -10T,  crack 
dimensions  2x2  cm  (after  Petrenko 
1992a). 
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Figure  25.  Schematic  of  experimental  configuration  used  in 
measurements  of  electric  fields  on  the  ice  surface  around  a 
stable  crack  (after  Petrenko  1992a). 


bandsaw  (1  mm  thick)  were  located  in  the  middle  of 
the  samples,  perpendicular  to  their  long  side.  Sample 
dimensions  were  15x6x6  cm^  and  the  crack  length 
varied  from  5  to  30  mm.  An  oscillating  sine  wave,  of 
frequency  from  1  to  300  Hz,  was  applied  to  the  sam¬ 
ples  by  an  MTS  testing  machine.  In  this  set  of  meas¬ 
urements,  a  sine  potential  difference  was  measured  by 
a  lock-in  amplifier  between  one  remote  electrode 
mounted  into  the  ice  and  one  electrode  placed  near  the 


Figure  26  shows  the  distribution  of  the  amplitude  of 
an  oscillating  electrical  potential  (p  measured  around 
an  artificial  mode  I  tearing  crack.  To  compare  these 
experimental  results  with  the  theory,  (p  was  calculated 
using  eq  51 


(53) 


crack  tip.  This  latter  electrode  in  some  cases  was 
mounted  into  the  ice  to  a  depth  of  2  mm,  by  careful 
heating  of  the  electrode,  or  was  attached  to  the  ice  sur¬ 
face  with  slight  pressure.  The  electrodes  were  made  of 
stainless  steel  wire  of  0.1 -mm  diameter. 


where  must  be  taken  from  eq  47  and  Oj  =  a4  =  6.2 
eV  (Evtushenko  et  al.  1987),  There  is  satisfactory 
quantitative  and  qualitative  agreement  between  theory 
and  experiment  in  the  vicinity  of  the  crack  tip  for  r  = 
1-10  mm.  At  r  =  1  mm,  (p  =  0.8  mV.  Yet,  at  larger  dis¬ 


tances  (p  remains  almost  constant  instead  of  decreasing 
to  0  as  the  theory  predicts  (9  r  r  is  the  distance 


Sample  dimensions  13  x  6.5  x  6.5  cm 
Crack  length  /  =  3  cm 
Width  d=  1  mm 
r=-10'’C 

Stress  amplitude  0.12  MPa 
Frequency  90  Hz 


Figure  26.  Amplitude  of  electric  potential  (pon  the  surface  of 
the  ice  sample  measured  in  the  configuration  shown  in  Fig¬ 
ure  25;  large  grains  of  ice  grown  from  deionized  water  (after 
Petrenko  1992a). 


to  the  tip).  This  discrepancy  could  arise  from  small  os¬ 
cillations  of  the  crack  surfaces  if  there  is  a  constant 
potential  difference  between  them. 

Similar  EME  from  cracks  were  observed  by  Thiel 
(1992).  He  recorded  EME  in  the  vicinity  of  ice  cracks 
created  by  breaking  rods  of  ice,  scouring  the  ice  sur¬ 
face,  applying  uniaxial  pressure  to  a  cylindrical  ice 
core  and  using  a  bore-hole  jack  in  an  ice  sheet.  The 
observations  were  made  in  the  frequency  band  from 
50  Hz  to  15  kHz,  using  a  two-channel  audio  cassette 
recorder  and  wire  probes  inside  or  outside  the  ice. 
Thiel  found  that  EME  events  correlate  with  acoustic 
emissions  events  that  arise  during  crack  formation. 

Field  experiments 

At  present  we  have  several  reports  on  observations 
of  EME  from  lake  ice,  glaciers  and  sea  ice  sheets.  Ka- 
churin  et  al.  (1979)  reported  the  detection  of  electrical 
and  magnetic  signals  from  ice  and  snow  under  field 
conditions.  Radio  waves  emitted  by  the  motion  of 
snow  avalanches  were  recorded  by  a  set  of  radio-fre¬ 
quency  receivers  working  in  the  frequency  ranges  of 
0. 1  to  30  Hz,  900  Hz  to  2  kHz,  2.5  MHz,  40  MHz  and 
760  MHz.  The  electrical  component  of  EME  was  re¬ 
corded  at  760  MHz,  and  the  magnetic  component  at 
the  other  frequencies.  The  sensitivity  of  the  detectors 
over  all  frequencies  was  1  to  2  |iV/m.  Broadband  EME 
with  a  maximum  near  1  kHz  were  recorded  for  ava¬ 
lanches  of  freshly  fallen  snow  and  near  2.5  MHz  for 
packed  snow.  The  amplitudes  of  the  electromagnetic 
signals  investigated  were  one  to  two  orders  of  magni¬ 
tude  higher  than  that  of  background  noise  in  the  same 
corresponding  range.  The  authors  did  not  provide  in¬ 


formation  either  on  the  distance  from  the  receiver  to  the 
avalanches  nor  on  the  types  of  antenna  used.  While  all 
the  mechanisms  of  EME  from  ice  considered  above  are 
applicable  for  numerous  cracks  in  a  snow  cover  during 
an  avalanche,  there  may  be  some  additional  phenome¬ 
na,  such  as  frictional  electrification,  responsible  for  the 
EME  from  snow  avalanches. 

The  same  authors  reported  detection  of  magnetic  sig¬ 
nals  captured  from  a  mountain  glacier  (Malyy  Azau  gla¬ 
cier  in  the  Caucasus  Mountains)  during  various  natural 
and  artificially  induced  dynamic  processes.  The  EME 
and  seismic  signals  were  recorded  in  the  0.1-  to  30-Hz 
frequency  range.  The  magnetic  component  was  cap¬ 
tured  with  a  single-turn  loop  with  a  diameter  of  100  m 
and  with  an  induction  coil  with  a  Permalloy  core.  Figure 
27,  taken  from  the  paper  by  Kachurin  et  al.  (1979), 
shows  the  correlation  between  EME  and  seismic  signals 
captured  from  the  glacier.  Later,  Kachurin  et  al.  (1984) 
reported  EME  from  freshwater  lake  ice  and  sea  ice 
sheets.  To  measure  the  horizontal  component  of  the 
magnetic  field,  they  used  an  induction  coil  with  a  reso¬ 
nance  frequency  of  10  Hz.  The  EME  from  lake  ice  were 
detected  during  air  cooling  at  night  and  during  ice  hum- 
mocking.  In  the  same  study  Kachurin  and  his  co-authors 
measured  quasi-static  (10-2  ^  1  Hz)  electrical  fields 
generated  in  first-year  sea  ice  during  its  loading.  When 
the  load  was  applied  they  detected  an  electrical  potential 
difference  of  a  few  millivolts  between  two  electrodes 
mounted  into  the  ice  at  a  distance  of  50  m  from  each 
other.  The  corresponding  bend  of  the  ice  sheet  was 
about  10"^  radian. 

They  also  made  a  very  interesting  attempt  to  detect 
sea  ice  EME  from  an  airplane  flying  at  an  altitude  of  100 
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Figure  27.  Synchronous  records  of  the  electromagnetic  (top)  and  seis¬ 
mic  (bottom)  signals  captured  from  natural  and  artificially  induced  gla¬ 
cier  tremors  on  17  August  1977  (after  Kachurin  et  al.  1979).  The  pass 
band  of  the  electromagnetic  and  seismic  signals  is  0.1  to  30  Hz. 
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1.  Fast  1-year  ice,  100-130  cm,  snow  20  cm, 
hummocks 

2.  Marine  ice  60-80  cm,  hummocks,  snow  20  cm, 
compression 

3.  Ice  fields  with  open  pools 

4.  Fresh,  snow-free  ice  30-40  cm,  compression 

5.  Large  open  pools,  ice-free  water,  tow  vessel 
released  and  180°  turn  made 

6.  Fresh  ice  with  open  pools 

7.  Ice  fields  with  open  pools 

8.  Passage  into  solid  ice. 


Figure  28.  Synchronous  3-day  en  route  measurements  of  electric  (top)  and  magnetic  (bottom)  fields  (after  Kachurin 
etal  1988). 


m  using  a  magnetic  ferrite  antenna  connected  to  a  selec¬ 
tive  amplifier.  The  amplifier  was  tuned  to  the  frequency 
of  100  kHz.  An  equivalent  sensitivity  of  the  measuring 
system  equaled  2  |iV/(mV/m),  and  an  equivalent  noise 
level  was  0.3  mV/m.  The  plane  flew  above  various  re¬ 
gions  such  as  open  water,  thin  ice  and  thick  ice  during 
different  wind  conditions.  Strong  wind  loads  the  ice, 
causing  intensive  cracking,  hummocking  and  deforma¬ 
tion.  They  found  that  the  strength  of  the  electromagnet¬ 
ic  field  increased  to  5-8  mV/m  above  the  regions  with 
intensive  ice  deformation,  while  without  wind  loading 
or  above  open  water,  variations  of  EME  did  not  exceed 
0.5  mV/m  as  measured  on  different  days. 

During  the  navigation  period  of  1986,  the  electrical 
and  electromagnetic  fields  arising  from  a  fracturing  of 
the  ice  cover  during  ship  passage  were  measured  from 
the  nuclear  icebreaker  A rA:n/:a  (Kachurin  et  al.  1988).  A 
dynamic  electrometer  was  used  for  the  measurements 
of  electrostatic  fields  in  the  range  of  10  to  10^  V/m.  The 
electrometer  was  on  a  movable  arm  mounted  11m 
above  the  waterline  at  the  bow  of  the  icebreaker.  Elec¬ 
tromagnetic  pulses  were  recorded  with  an  induction 
transducer  having  a  ferrite  core  with  a  resonance  fre¬ 
quency  of  about  1  kHz.  The  signals  were  amplified  by  a 
selective  microvoltmeter  in  the  frequency  range  from 


200  Hz  to  2  kHz.  The  researchers  picked  this  frequency 
range  to  minimize  the  level  of  natural  and  industrial 
electromagnetic  interference  and  to  allow  subsequent 
satellite  measurements  through  the  ionospheric  radio 
propagation  window  in  the  ELF  range.  An  equivalent 
ambient  noise  was  about  10“^®  T  while  the  signals  usu¬ 
ally  measured  from  ice  fracture  exceeded  10“^  T. 

Figure  28  shows  the  results  of  synchronous  measure¬ 
ments  of  the  electrical  and  electromagnetic  fields  made 
in  ice  sheets  of  various  thicknesses  and  ages.  Icebreaker 
velocity  and  electrical  and  magnetic  field  strengths  cor¬ 
related  over  a  much  shorter  time  scale  of  several  min¬ 
utes.  The  electrical  and  magnetic  fields  measured  in 
these  experiments  may  be  caused  by  both  EME  from 
cracks  and  frictional  electrification  of  ice  and  snow  gen¬ 
erated  by  the  icebreaker  passing  through  the  snow-cov¬ 
ered  ice. 

Petrenko  and  Gluschenkov  (1995)  reported  field 
measurements  of  EME  from  individual  cracks  in  fresh¬ 
water  lake  ice,  and  Gluschenkov  and  Petrenko  (1993) 
did  the  same  with  first-year  sea  ice.  The  electrical  com¬ 
ponent  of  the  EME  from  cracks  was  measured  with  an 
experimental  setup  shown  schematically  in  Figure  29. 
A  10-m-long  dipole  antenna  was  stretched  horizontally 
1  m  above  the  ice  surface.  The  antenna  was  connected  to 
an  input  circuit  with  a  high  input  impedance  (10^^  Q) 
preamplifier.  The  high  input  impedance  together  with 
the  electrical  capacitance  of  the  antenna  ('-lO"^^  F)  en¬ 
abled  us  to  measure  electrical  fields  in  a  wide  frequency 


Figure  29.  Experimental  configuration  used  in  meas¬ 
urements  of  EME  from  cracks  in  lake  and  sea  ice  (after 
Gluschenkov  and  Petrenko  1993). 
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Figure  30.  Typical  signal  captured  from  a  thermo^ 
crack  grown  in  lake  ice  (after  Gluschenkov  and 
Petrenko  1993).  T 


range,  from  about  10“^  Hz  to  1  MHz.  Additional  filter¬ 
ing  was  used  to  cut  off  the  common  industrial  (60  Hz) 
and  ionospheric  (/  >  100  kHz)  noise. 

Most  of  the  signals  captured  from  identified  cracks 
were  very  similar  to  those  found  in  the  laboratory  (Pe¬ 
trenko  1992a,  Fifolt  et  al.  1993).  A  typical  signal  cap¬ 
tured  from  a  thermocrack  in  freshwater  columnar  lake 
ice  is  depicted  in  Figure  30.  The  crack  appeared  10  m 
from  the  antenna  and  split  several  single-crystalline 
columns  of  ice.  As  in  the  laboratory  experiments,  the 
oscillation  of  the  electrical  field  can  be  interpreted  in 
terms  of  an  oscillating  electrical  charge  on  the  crack 
surfaces. 

An  electrical  signal  captured  by  a  dipole  antenna 
from  a  larger  crack  splitting  a  cold  (-35°C)  plate  of  sea 
ice  is  shown  in  Figure  31.  The  crack  surface  was  nor¬ 
mal  to  the  direction  of  the  gradient  of  salinity  of  the  ice, 
and  hence  the  crack  grew  perpendicular  to  the  frozen-in 
intrinsic  electrical  field.  Again,  as  in  corresponding  lab¬ 
oratory  experiments,  the  signal  increases  during  crack 
propagation  and  then  relaxes  after  crack  arrest. 


Figure  31.  Electromagnetic  emission  from  a  crack  in 
sea  ice  (after  Petrenko  and  Gluschenkov  1995).  The 
crack  was  60  cm  long  and  10  cm  wide.  Area  1  corre¬ 
sponds  to  the  crack  propagation  and  area  11  corre¬ 
sponds  to  electrical  relaxation  of  the  electric  field  after 
crack  arrest. 

ELECTROPLASTIC  EFFECTS  IN  ICE 

Several  electromechanical  phenomena  in  ice  are 
associated  with  the  motion  of  electrically  charged  dis¬ 
locations.  They  are  all  of  quite  small  magnitude.  Nev¬ 
ertheless,  their  significance  is  determined  by  their  con¬ 
tribution  to  the  study  of  the  physical  mechanisms  that 
govern  motion  of  dislocations  and,  hence,  plastic  de¬ 
formation  of  ice.  Electrically  charged  dislocations 
were  found  in  many  materials  such  as  ionic  crystals 
(see  review  by  Whitworth  1975);  covalent  semicon¬ 
ductors  (see  review  by  Alexander  and  Teichler  1991); 
and  crystals  with  mixed  covalent-ionic  bonding  (see 
review  by  Osip’yan  et  al.  1986).  In  these  materials  the 
electrical  charge  on  a  dislocation  core  appears  when  a 
dislocation  captures  either  charged  point  defects  (such 
as  vacancies  and  interstices  in  ionic  crystals)  or  elec¬ 
trons  (or  holes),  as  happens  in  semiconductors  or  when 
both  point  defects  and  electrons  (holes)  are  captured. 
Significant  theoretical  efforts  have  been  applied  to 
clarify  the  formation  of  the  dislocation  charges  in  ionic 
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Figure  32,  Section  in  the  (ITOO)  plane  of  a  60  dislocation 
on  a  plane  of  the  shuffle  set  in  the  structure  of  ice,  illustrat¬ 
ing  the  dangling  bond  in  the  core. 


materials  and  semiconductors.  Though  there  are  no 
theoretical  papers  on  the  electrical  charge  of  disloca¬ 
tions  in  ice,  some  reasons  for  formation  of  such  a 
charge  are  obvious. 

First,  the  charge  may  originate  from  electrical  ac¬ 
tivity  of  dangling  hydrogen  bonds  in  the  dislocation 
core.  Figure  32  illustrates  a  dangling  bond  in  the  core 
of  a  60°  dislocation.  A  dangling  bond  may  or  may  not 
have  a  proton  on  it.  In  the  first  approximation,  the 
probability  of  finding  a  proton  on  a  dangling  bond  is 
50%,  since  when  a  regular  hydrogen  bond  breaks  into 
two  dangling  bonds,  one  proton  is  shared  between  two 
bonds.  It  is  generally  acknowledged  that  a  regular  hy¬ 
drogen  bond  is  electrically  neutral.  The  dangling 
bonds  of  a  dislocation  core  introduce  electron  energy 
levels  or  narrow  energy  bonds  located  in  the  forbidden 
band  of  ice.  Depending  on  the  relative  position  of  a 
Fermi  level  and  the  dislocation  levels  (bands),  the  ex¬ 
change  of  electrons  between  the  dangling  bonds  and 
the  bulk  ice  will  result  in  a  net  negative  or  positive 
charge  of  the  dislocation.  It  is  not  known  at  present  if 
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Figure  33.  Directions  of  the  principal  stresses  in  the 
vicinity  of  an  edge  dislocation.  The  dislocation  line  is 
directed  perpendicularly  to  the  figure*  s  plane. 


such  dangling  bonds  may  exchange  protons  with  an 
ice  bulk. 

Second,  a  mechanism  that  may  result  in  a  net  dislo¬ 
cation  charge  is  strong  elastic  interaction  between  a 
protonic  charge  carrier — i.e.,  ions  and  Bjerrum  de¬ 
fects — and  strain  generated  by  a  dislocation.  Figure  33 
illustrates  stress  patterns  in  the  vicinity  of  an  edge 
dislocation.  As  is  seen  from  this  figure,  the  dislocation 
generates  compression  in  the  upper  quadrants  and  ten¬ 
sion  in  the  lower  ones. 

The  specific  volumes  of  protonic  charge  carriers 
(ions  and  Bjerrum  defects)  exceed  that  of  a  water  mol¬ 
ecule  (Evtuschenko  et  al.  1987,  Evtuschenko  and 
Petrenko  1991).  Thus,  all  these  charged  defects  should 
accumulate  in  the  lower  quadrant  and  escape  from  the 
upper  one.  Figure  34  shows  schematically  the  change 
in  defect  activation  energy  as  a  function  of  the  a:  coor¬ 
dinate,  also  shown  in  Figure  33.  Since  for  the  protonic 
defects  in  ice  aj  ranges  from  1.2  to  6.8  eV  and  Sa  is 
about  0.1  in  the  vicinity  of  a  dislocation  core,  the  bind¬ 
ing  energy,  oq  Eij,  of  the  defects  on  an  edge  dislocation 


Figure  34.  Change  in  activation  energy  aza  of  a  charge 
carrier  in  ice  in  the  vicinity  of  an  edge  dislocation. 


is  about  tens  of  electron  volts  and  is  comparable  with  the 
defects’  activation  energy.  That  should  result  in  signifi¬ 
cant  buildup  of  defects  near  a  dislocation  and  in  total 
electrical  charge  if  a;  differs  for  positive  and  negative 
defects.  Such  a  mechanism  of  dislocation  charge  cre¬ 
ation  was  first  suggested  and  estimated  in  Petrenko  and 
Ryzhkin  (1986a, b),  but  has  never  been  accurately  calcu¬ 
lated. 

Electrically  charged  dislocations  can  manifest  them¬ 
selves  in  three  different  electromechanical  phenomena. 
First,  their  motion  during  plastic  deformation  can  result 
in  an  electrical  current,  a  so-called  dislocation  current. 
Second,  application  of  an  external  electrical  field  may 
cause  the  dislocation  motion  and,  hence,  small  plastic 
deformation.  Third,  the  introduction  of  dislocations 
may  change  concentrations  of  charge  carriers  and,  as 
a  result,  cause  a  change  in  electrical  conductivity  and 
dielectric  permittivity  of  ice.  We  will  consider  below  all 
three  groups  of  phenomena. 

Dislocation  currents  in  ice 

Dislocation  currents  in  ice  were  measured  and  used 
for  determination  of  the  linear  density  of  dislocation 
charge  q  (Petrenko  and  Whitworth  1983).  We  used 
specimens  of  very  pure  single  crystals  of  ice  in  which 
the  majority  of  edge  dislocations  of  one  mechanical 
sign  was  introduced  by  preliminary  bending,  as  shown 


in  Figure  35a.  A  smaller  ice  sample  was  then  cut  off 
the  bent  specimen  for  tensile  tests  (Fig.  35b  and  c).  The 
electrodes  were  formed  from  an  Hg-In  amalgam, 
which  adheres  well  to  the  surface  of  ice  but  remains 
liquid  down  to  -30°C,  so  that  it  does  not  interfere  with 
deformation.  Dislocation  currents  were  measured  with 
an  electrometer  at  7  =  -20°C  during  tensile  creep  of 
the  ice  samples.  Atypical  recording  of  both  the  current 
/  and  the  tensile  elongation  Af  are  shown  in  Figure 
35d.  The  abrupt  changes  in  the  current  A/  were  caused 
by  the  starts  and  stops  of  charged  dislocations.  It  was 
shown  that  q  can  be  calculated  as 


where  hz=  component  of  the  Burgers’  vector  along 
the  length  of  the  specimen 
d=  interelectrode  distance 
AI=  dislocation  current 
f  p  =  rate  of  the  tensile  plastic  deformation. 

The  dislocation  charge  may  be  expressed  in  terms  of 
qale,  i.e.,  in  the  number  of  elementary  charges  per  in- 
termolecular  distance  a.  We  found  that  the  dislocation 
charge  in  ice  was  positive,  with  an  absolute  magnitude 
of  qaJe  =  2  x  10“3.  That  corresponded  to  two  proton 


a.  Bent  crystal  containing  an  excess  of  disloca¬ 
tions  of  one  sign.  When  deformed  by  tension  the 
majority  of  dislocations  move  toward  the  upper 
electrode  Ej. 


b.  How  the  tensile  specimen  is  cut  from  the  bent  crystal. 


c.  Ice  specimen  mounted  for  ten¬ 
sile  deformation  experiment. 
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d.  Recordings  of  change  of  length  A£  and 
current  I  during  tensile  deformation.  Load  is 
applied  and  removed  at  times  indicated  by  the 
arrows.  A£  was  measured  in  micrometers. 


Figure  35.  Dislocation  currents  in  ice  (after  Petrenko  and  Whitworth  1983). 
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charges  per  thousand  intermolecular  distances  along 
the  dislocation  core.  This  charge  magnitude  was  esti¬ 
mated  as  a  lower  limit  because  of  the  use  of  blocking 
(non-ohmic)  electrodes,  so  that  the  charge  flow  within 
the  ice  produced  polarization  at  the  surface  and  the  full 
current  may  not  be  recorded  by  the  electrometer. 


Motion  of  charged  dislocations 
in  an  electrical  Held 

Itagaki  (1970)  used  a  different  technique  to  deter¬ 
mine  the  dislocation  charge  density  q.  Using  X-ray  to¬ 
pography  to  observe  dislocation  lines  in  ice,  he  found 
that  some  dislocations  vibrated  when  an  external  ac 
electrical  field  was  applied.  At  the  low-frequency  limit, 
the  amplitude  of  vibrations  Tj  depends  only  on  the 
length  of  a  dislocation  segment  external  electrical 
field  E  and  the  tension  of  a  dislocation  line  T^.  Namely 


ri  = 


87^ 


(55) 


This  equation  enables  us  to  estimate  qif  and  E  are 
known. 

Itagaki  measured  the  amplitudes  of  vibrations  T)  in 
pure  single  crystals  of  ice  taken  from  the  Mendenhall 
glacier.  The  ac  field  strength  ranged  from  3  to  600 
V/cm,  and  three  frequencies  of  60,  1.33  and  0.33  Hz 
were  used.  The  linear  density  of  a  dislocation  charge 
estimated  by  Itagaki  ranged  from  qa/e  =  3.6  X  10“3  to 
qa/e  =  2  x  lO”^. 

Later,  Joncich  et  al.  (1978)  tried  to  explore  a  similar 
idea  to  estimate  q  for  dislocations.  They  applied  a 
strong  electrical  field  across  a  low-angle  title  boundary 
consisting  of  edge  dislocations.  Their  tests  were  per¬ 
formed  at  r  =  -15°C.  No  measurable  displacement  of 
the  boundary  was  observed.  Based  on  their  experi¬ 
mental  precision,  Joncich  et  al.  concluded  that  qa/e 
cannot  exceed  3.3  x  ICL^,  Since  they  used  a  dc  electri¬ 
cal  field  and  blocking  electrode,  the  real  E  in  the  ice 
bulk  in  their  experiment  could  be  lower  than  they  ex- 
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pected.  This  should  increase  the  estimated  upper  limit 
of  qa/e  somewhat. 

We  can  conclude  that  there  is  no  significant  discrep¬ 
ancy  among  different  authors  and  the  techniques 
applied  and  that  most  probably  the  magnitude  of  qa/e 
should  be  somewhere  in  the  range  from  2  x  1()~3  to  2  x 
10“2.  If  the  above-described  accumulation  of  point  de¬ 
fects  in  the  region  of  tension  near  a  dislocation  core  is 
responsible  for  the  dislocation  charge,  then  the  small 
linear  density  of  this  charge  does  not  reflect  the  real 
number  of  the  defects  captured,  but  rather  shows  a  dif¬ 
ference  in  concentrations  of  positive  and  negative 
defects  near  the  dislocation. 

Effect  of  static  electrical  field  on  ice  creep 

In  ionic  crystals  and  some  semiconductors,  qa/e  is 
so  large  (0.1— 1.0)  that  the  application  of  high  external 
electrical  fields  of  £  =  10^  V/cm  is  able  to  accelerate  or 
halt  creep,  depending  on  the  relative  directions  of  qE 
and  the  predominant  velocity  of  dislocations  (Whit¬ 
worth  1975,  Osip’yan  et  al.  1986).  Such  an  effect  of 
electrical  fields  on  plastic  deformation  in  ice  has  not 
been  observed  so  far.  However,  the  application  of  even 
smaller  static  electrical  fields  of  10^  V/cm  can  signifi¬ 
cantly  suppress  plastic  creep  of  thin  ice  specimens  hav¬ 
ing  a  thickness  of  1  mm  or  less  (Petrenko  and  Schulson 
1993). 

We  examined  the  effect  of  dc  and  ac  electrical  fields 
on  the  plastic  deformation  of  pure  and  doped  (with  HF 
and  KOH)  single  crystals  of  ice  in  the  temperature 
range  of  -10  to  ^5°C.  The  basal  plane  (0001)  was  a 
slip  plane  for  dislocations.  We  found  that  a  relatively 
small  (about  10—50  V)  dc  voltage  applied  across  thin 
single  crystals  decreased  both  the  creep  strain  rates 
and  the  high-frequency  conductivity  of  ice  a„.  Figure 
36  depicts  the  dependence  of  the  creep  strain  rate  on  the 
dc  voltage  applied  to  a  50-|im  thick  ice  specimen.  At 
the  same  time,  the  application  of  an  ac  electrical  field  of 
the  same  strength  as  the  dc  field  or  the  application  of  a 
dc  field  to  thick  specimens  (1-3  mm)  did  not  reveal 
changes  in  the  creep  rate.  We  felt  that  the  observed 
changes  in  the  plastic  deformation  were  ascribable  to  a 
decrease  in  the  concentrations  of  Bjerrum  defects 
caused  by  the  strong  static  electrical  fields.  The  con¬ 
centration  of  Bjerrum  defects  is  proportional  to  o„  and 
can  be  easily  monitored. 


Figure  36.  Dependence  of  the  relative  change  in  strain 
rate,  e  on  the  applied  dc  voltage  V  (after  Petrenko  and 
Schulson  1993).  i  is  the  strain  rate  a/  V  =  0.  r=  -20®C, 
sample  thickness  L  =  50  \xm. 
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Effect  of  plastic  deformation  on 
electrical  properties  of  ice 

Dislocations  introduced  during  the  process  of  plas¬ 
tic  deformation  can  affect  the  electrical  properties  of 
ice  (as  well  as  other  materials)  in  a  variety  of  ways. 
They  can  generate  new  point  defects  such  as  vacan¬ 
cies,  interstitials  and  Bjerrum  defects;  they  can  also 
capture  point  defects  by  their  elastic  field  and  make 
defects  immobile.  All  this  can  change  the  charge  carri¬ 
er  concentrations.  But  this  is  not  the  whole  story  yet. 
Owing  to  their  long-range  elastic  and  electrical  fields, 
dislocations  can  effectively  scatter  charge  carriers.  As 
a  result  the  charge  carrier  mobility  usually  decreases 
after  plastic  deformation.  Moreover,  charged  disloca¬ 
tions  may,  in  principle,  move  under  the  action  of  these 
electrical  fields,  contributing  to  both  electrical  conduc¬ 
tivity  and  dielectric  permittivity. 

The  first  known,  but  unsuccessful,  attempt  to  find 
the  effect  of  plastic  deformation  on  the  electrical  prop¬ 
erties  of  ice  was  made  by  Brill  and  Camp  (1957).  They 
reported  no  changes  induced  by  plastic  deformation  in 
the  dielectric  properties  of  ice  at  1  kHz.  On  the  con¬ 
trary,  Higashi  (1969)  reported  preliminary  results  ob¬ 
tained  by  Mae,  who  found  a  very  profound  increase  in 
low-frequency  conductivity  and  dielectric  permittivity 


of  ice  that  had  undergone  15  and  30%  plastic  deforma¬ 
tion  (Fig.  37).  These  changes  were  confirmed  later  by 
Mae  and  Higashi  (1973),  after  they  conducted  very  de¬ 
tailed  and  careful  investigations  on  the  effect  of  plastic 
deformation  on  ice  electrical  conductivity  and  dielectric 
permittivity  in  wide  temperature  (-10  to  -100°C)  and 
frequency  (0.5  to  10^  Hz)  ranges.  During  plastic  defor¬ 
mation  of  pure  single  crystals  of  ice,  the  dislocation 
density  in  their  experiments  rose  from  10'^  to  lO"^  cm"^. 
Mae  and  Higashi  found  that  plastic  deformation  left  un¬ 
changed  the  high-frequency  conductivity  ctco  but  in¬ 
duced  significant  changes  in  low-frequency  conductivi¬ 
ty  Og.  This  can  be  easily  explained,  as  the  plastic  defor¬ 
mation  generated  a  concentration  of  new  charge  carriers 
that  was  comparable  with  the  minority  charge  carrier 
concentration  but  was  much  less  than  the  majority 
charge  carrier  concentration.  That  would  also  explain 
the  failure  of  Brill  and  Camp,  who  measured  E  in  the 
region  of  majority  charge  carrier  fi-equency. 

A  typical  change  in  Og  found  by  Mae  and  Higashi, 
after  introducing  a  dislocation  density  of  lO'^/cm^,  was  2 
X  10-*0/Q  -  cm.  Notice  that  this  number,  when  recalcu¬ 
lated  into  the  number  of  elementary  charges  per  inter- 
molecular  length  on  a  dislocation,  gives  5  x  10“3,  which 
is  in  the  range  of  the  experimental  magnitude  of  qa/e. 
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Figure  37.  Effect  of  compressive  strain  on  the  electric  properties 
of  ice  single  crystals,  measured  at  1  Hz  (after  Higashi  1969). 
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Mae  and  Higashi  also  analyzed  separately  changes 
in  the  mobility  and  concentration  of  minority  charge 
carriers  caused  by  plastic  deformation.  They  used  for 
this  purpose  their  experimentally  obtained  frequency 
dependencies  of  a  and  e,  but  exploited  a  theory 
by  MacDonald  (1953),  which  is  not  applicable  to  ice. 
Their  results  showed  unrealistic  magnitudes  for  an 
ion  mobility  of  0.6  cm^A^  s.  Nevertheless,  when  the 
same  experimental  technique  is  combined  with  a  spe¬ 
cific  theory  for  ice  (Petrenko  and  Ryzhkin  1984a), 
it  gives  a  correct  value  for  mobility  of  ions,  of  a  few 
10“^  cm^A^  s  (Zaretskii  1991). 

Noll  (1978)  studied  the  influence  of  plastic  defor¬ 
mation  on  electrical  properties  during  and  after  defor¬ 
mation.  He  measured  o  and  e  at  7  -  -3  and  -10°C  in  a 
frequency  range  from  10~2  to  10^  Hz,  a  strain  rate  range 
from  2x  10“^  to  10“^  s“^  and  a  strain  range  from  0.03 
to  0.6.  He  found  that  plastic  deformation  reduces  the 
electrical  conductivity  in  the  space-charge  dispersion 
frequency  range  (Fig.  38).  The  decrease  depended 
strongly  on  the  strain  rate,  while  the  strain  magnitude 


Figure  38.  Effect  of  plastic  deformation  on  dielectric  per¬ 
mittivity  t  and  conductivity  a  (after  Noll  1978).  Solid  lines 
were  obtained  before  deformation,  dashed  lines  were  obtained 
during  plastic  deformation  at  £p  =  3.2  x  IQ-^  ^  and  £p  =  0.274. 
T=  210.2  K. 


had  little  effect.  After  deformation  takes  place,  recov¬ 
ery  processes  have  been  observed  that  restore  the  orig¬ 
inal  conductivity  if  the  total  strain  was  small.  Rapid 
and  large  deformation,  however,  causes  permanent 
changes  in  the  conductivity  and  dielectric  permittivity. 

Itagaki  (1978,  1982,  1983)  conducted  a  very  inter¬ 
esting  series  of  experiments  and  theoretical  calcula¬ 
tions  aimed  at  determining  the  contribution  of  electri¬ 
cally  charged  dislocations  to  the  electrical  properties 
of  ice.  In  1978  he  carried  out  an  exciting  experiment  in 
which  he  compared  the  electrical  properties  of  disloca¬ 
tion-free  areas  of  ice  pure  microcrystals  with  areas 
having  a  high  density  of  dislocations.  He  used  hoar¬ 
frost  crystals  as  the  specimens  and  X-ray  topography 
to  visualize  dislocations  in  them.  Use  of  tiny  liquid- 
mercury  electrodes  prevented  damage  to  the  ice  crys¬ 
tals.  Itagaki  did  not  observe  the  Debye  dispersion  of  a 
and  £in  dislocation-free  regions,  while  it  was  present 
in  the  regions  with  high  dislocation  density  and  in  the 
initially  dislocation-free  regions  after  introduction  of 
dislocations  by  surface  scratching.  On  the  basis  of 
these  observations,  he  concluded  that  the  main  reasons 
known  for  the  polarization  process  in  ice,  the  De¬ 
bye  dispersion  and  hence  the  high-frequency  con¬ 
ductivity  Goo,  are  related  to  the  motion  of  charged 
dislocations  in  an  ac  electrical  field. 

He  later  developed  a  theoretical  model  support¬ 
ing  this  idea  (Itagaki  1982,  1983).  However,  Ita¬ 
gaki ’s  model  was  not  accepted  by  the  ice  physics 
community.  While  dislocations  in  ice  can  contrib¬ 
ute  to  and  modify  the  electrical  properties  of  ice, 
they  cannot  account  completely  for  those  electri¬ 
cal  properties.  There  may  be  other  reasons  for  the 
absence  of  Debye  dispersion  found  by  Itagaki  in 
very  thin  and  pure  specimens  of  ice.  One  possibil¬ 
ity  is  that  the  specimen  thickness  was  comparable 
with  or  less  than  a  screening  length.  Conventional 
theories  of  electrical  properties  of  ice  (Jaccard 
1964)  predict  an  absence  of  the  Debye  dispersion 
in  this  case  (Petrenko  and  Ryzhkin  1984a). 
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